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Stability Theary, Aeraaacuatioa , Shear layera 
Excited Waves In Shear Layers 

Summary 

The generation of instability waves in free shear layers is in- 
vestigated. The model assumes an infinitesimally thin shear layer 
shed from a semi-infinite plate which is exposed to sound exci- 
tation. For this model * it is shown that only forced instability 
waves and a Kutta condition at the trailing edge are possible. 

The shear layer excitation by a source further away from the 
plate edge in the downstream direction is very weak while up- 
stream from the plate edge the excitation is relatively effi- 
cient. A special solution is given for the source at the plate 
edge. Any type of source further away from the plate edge creates 
a parabolic pressure field near the edge. For this latter, fairly 
general case, a reference quantity is found for the magnitude of 
the excited instability waves. The theory is then extended to 
two streams on both sides of the shear layer having different 
velocities and densities. Furthermore, the excitation of a shear 
layer in a channel is calculated. The limitations to the theo- 
ry and some aspects related to experiments are discussed. In 
particular, for a comparison with measurements, numerical com- 
putations of the velocity field outside the shear layer have 
been carried out. 

Stabilittitatheorie , Aeroakiutik, Schereehiohten 
Angeregte Wellen in Scherschichten 

Ubersicht 


Die Erzeugung von InstabilitMtswellen in freien Scherschichten 
wird untersucht. Dabei wird eine unendlich dllnne Scherschicht 
betrachtet, die am Ende einer halbunendlichen Platte entsteht. 
Diese Scherschicht wird akustisch angeregt. FUr dieses Modell 
wird gezeigt, daB nur erzwungene InstabilitMtswellen mbglich 
sind und daB an der Plattenendkante die Kuttasche AbfluBbedin- 
gung gilt. Die Anregung der Scherschicht durch eine Quelle strom- 
ab der Plattenkante ist schwach, wShrend die AnregUng von strom- 
auf relativ wirksam ist. Filr die Anregung direkt an der Platten- 
kante wird eine spezielle Lbsung angegeben. Jede Art von Quellen 
in groBerer Entfernung von der Plattenkante erzeugt ein parabo- 
lisches Druckfeld in KantennShe. Zu diesem letzteren, recht all- 
gemeinen Fall wird eine ReferenzgrSBe fttr die Amplitude der In- 
stabilitatswellen angegeben. Die Theorie wird ausgeweitet auf 
zwei Stromungen auf beiden Seiten der Scherschicht mit verschie- 
denen Geschwindigkeiten und Dichten. AuBerdem wird die Anregung 
einer Scherschicht in einem Kanal berechnet. Die Grenzen der 
Theorie sowie einige experimentelle Aspekte werden diskutiert. 
Zum Vergleich mit Messungen wird das Geschwindigkeitsfeld auBer- 
halb der Scherschicht numerisch berechnet. 
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Let us first consider a picture of the interaction process which 
we are going to investigate. In Fig. 1 a photograph of the evo- 
lution of a plane free shear layer can be seen. The shear layer 
is forced by acoustic excitation into wave motion and vortex 
roll-up. The lower half of the picture shows a flow from left to 
right whereas the upper half shows a very slow entrainment flow 
in the same direction. A smoke filament is injected into the 
shear layer from the low-Bpeed side. A splitter plate separating 
the two streams ends at the left hand side of the picture. Around 
the end of this splitter plate a fluctuating pressure field is 
created by vibrating plates outside the flow. 


splitter plate 

smoke injection 


entrainment flow 




|=| 1cm 


-► mean flow 


Fig. 1 Free shear layer evolution with acoustic forcing. For 
more convenient reproduction, the colours have been inverted so 
that the white sulfuric acid smoke filament is shown in black. 
Mean flow velocity 6.4 m/s, entrainment velocity « 0.6 m/s, fre- 
quency 63 Hz, forcing reference velocity Ap 12 /pVw01 = 0.354 m/s 
RMS (definition see section 2.1.), Strouhal number S Q * 0.0040. 

Although impressive, pictures like Fig. 1 have been observed many 
times in the past. In particular, the observation of the sensi- 
tivity of flames to sound has been made by Leconte [1] as early 
as 1858. Tyndall [2] has shown in 1867, that also jets without 
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combustion are sensitive to sound. Experimental observations like 
these were already known to Lord Kelvin , von Helmholtz and Lord 
Rayleigh when they laid the foundations to the stability theory 
of flows. Lord Rayleigh , e. g, , refers to sound-sen3itive smoke 
jets in his 1880 paper “On the stability, or instability, of 
certain fluid motions" t3). Since that time, stability theory 
has used a model assuming a parallel flow extended from x * -« 
to x » +« which is disturbed by a wave motion with constant 
magnitude in streamwise direction and growing with time. This 
approach is perfectly valid to answer the question whether or not 
a flow is unstable. However, it is not appropriate to use this 
model if an instability wave is driven by an oscillatory excita- 
tion which does not vary with time. In 1962, it has be^n sug- 
gested by Gaster [4] and Watson [5J , that in this latter case a 
stability analysis with spatially growing waves is the only 
reasonable approach. For free shear layers, the spatial stabili- 
ty model has been verified by Michalke [6] and Freymuth [7], In 
addition, it has become evident from Freymuth 's experiments, that 
the magnitude of the spatial instability waves is linearly de- 
pendent on the exterior forcing by, e.g., an acoustic field. 

In spite of the numerous experiments with artificially forced 
shear layers and jets which have been carried out in the mean- 
time, the mechanism describing how the perturbations were intro- 
duced into the flow was not yet understood.*) About a decade 
ago, however, a few theoretical papers appeared which included 
a semi-infinite plate into the stability analysis of a thin free 
shear layer. The introduction of such a solid surface into the 


*) Also some theoretical papers on this issue have inherent 
problems, e.g., Tam's paper (8) discusses the excitation of a 
shear layer by an acoustic beam. His model assumes a beam width 
of typically 20 t 60 times the momentum thickness of the shear 
layer. On the other hand, the range of unstable waves in the 
shear layer requires a relatively low excitation frequency. This 
would cause the beam width to be about 1/100 of the acoustic 
wavelength. Unfortunately, it is not possible to form an acous- 
tic beam narrower than several acoustical wavelengths. 
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analysis turned oat to bo of crucial importance in the under* 
standing of the coupling between sound waves and instability 
waves, 'i'lio first of those papers [91 , albeit mathematically cor- 
rect. led to some confusion of the physics. Later, a few other 
theoretical papers were published {10, 11} which utilised (as 
{9}) a Fourier transform procedure (i.e,, the Wiener- liopf tech* 
niquo) , whereas a simpler approach using symme try /a symmetry 
arguments {12, 13) turned out to bo more olognnt but also less 
versatile in its applicability, in the present investigation we 
will extend and exploit this latter approach l 12) so that theo- 
retical data are available which can bo compared to experiments. 
Any testing of the new theories {10 * 13) had not boon carried 
out previously because the previous theories were presented in a 
way which could not bo verified easily by experimentalists, 
indeed, wo think that it is dosirablo in this field to track 
down by individual experiments ail basic assumptions and findings 
of the first simple theories before going on to more complex 
problems. The rewarding aspect of this ’’old fashioned" type of 
research is, that simple analytic solutions are obtained which 
can be checked with (moderately) simple experiments. 

There are different motivations for the present investigations 

(i) prom the viewpoint of scientific curiosity it is interest- 
ing how perturbations are introduced into a shear flow. 
These perturbations cause instability waves which, in turn, 
load to turbulence. This development can be seen in Fig. 1, 
The interaction region is found close to the lip of the 
splitter plate. The emerging instability waves are ampli- 
fied rapidly in the downstream direction. As usual in flow 
visualisation pictures, the linear growth regime cannot be 
seen clearly, because the displacements of the shear layer 
are still too small to be visible, Since the amplification 
rate of the instability waves is rather high (about 500 
times within the first 10 cm of Fig. 1), the onset of non- 
linearities characterized by the vortex roll-up occurs 
rather suddenly (see also (6, 7) ) . 
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(ii) A behaviour very similar to instability wave evolution 

(orderly structures) lias been found also in turbulent jets 
[14, 15] export to acoustic excitation. On the other hand, 
similar orderly structures have been found also in "un^ 
excited jet3« [16], Recently, however, it has been shown, 
that at least in turbulent shear layers, orderly structures 
become highly organized and dominant if an involuntarily 
generated, mostly unnoticed and sometimes inaudible acoustic 
excitation does exist [17, 18), 

(Hi) 'l'he investigations on acoustically excited turbulent jets 
have led to findings of practical importance, it has been 
found, that broad band jet noise can be amplified by a pure 
tone acoustic excitation from inside the nozzle [19 * 21), 

Because of its relevance to aircraft noise generation this 
observation has triggered a number of further investiga- 
tions [22 t 24). it has been recognized, that tho acoustic 
excitation of stability waves plays a crucial rfile in this 
mechanism (19, 21, 24], 

(iv) Tho change in acoustic radiation of turbulent jets duo to 
acoustical excitation is accompanied by quite substantial 
changes in the flow pattern. An amplification of tho jet 
noise corresponds to an enhanced mixing and an increased 
pertubation level in tho near field of the jet [21,25] , 

On the other hand, it has boon shown, that the spreading 
rate of a jet can also be reduced by acoustic control of 
the shear layer [25, 20). 

(v) whereas the above-mentioned broad band jet noise amplifi- 
cation is basically a nonlinear effect, another related 
interaction has been found. It has been shown, that the 
generation of instability waves at an edge can extract 
energy from a sound field [27 1 30). This absorption 
effect is (at least for small perturbations of a mean flow) 
a linear effect. It can lead to quite dramatic sound ab- 
sorptions, in particular at low frequencies. To give an 
example; Consider the sound transmission through a nozzle 
and a jet ah, say, a jet Mach number of W * 0.3 and a 
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dimensionles frequency ~ * 0.1, For these parameters*) 

“O 

only 1 % of the sound power transmitted through the nozzle 
exit will be transmitted into the far field! 

This paper, however, will not provide ary detailed prediction on 
any of the three latter effects. Nevertheless, the basic 
underlying interaction mechanism between sound waves and free 
shear layers will be dealt with and we will be able to predict, 
e.g., the perturbation magnitude introduced into the shear layer 
shown in Fig. 1. In addition, a well-defined reference quantity 
representing the acoustic excitation will be given sv that future 
experimental investigators will be able to replace arbitrarily 
taken reference pressures, velocities, or even loudspeaker input 
voltages by a more general and appropriate reference quantity. 

This paper will only deal with theoretical aspects of the shear 
layer excitation. There is evidence from our recent experiments 
that the present theory is valid in the expected range. 

2. One stream model 


Fig . 2 shows the simplified configuration which will be 
modelled mathematically. The geometry is basical . >..he same as 
in Fig. 1. We will, however, assume that there is no ^low 
above the shear layer (this condition will be released in 
section 5) . The acoustic field is assumed to be created by a 
pulsating source outside the shear layer in the fluid at rest. 

The following simplifying assumptions are introduced: 

1 . plane flow 

2. parallel mean flow 

3. linearized problem 

4. all fluctuating quantities harmonic in time , i . e . , “e"^ 1 

5. ioviscid flow 

6. incompressible flow 

7. infinitesimally thin shear layer 

*) R is the radius of the nozzle, oi is 2 tt times the frequency 
and a Q is the sound speed. 
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Fig. 2 Configuration of the analytic model, 


The first five simplifications are common in the stability theory 
of free jets. The assumption of an inviscid flow works quite well 
at sufficiently high Reynolds number [6, 7) . The sixth assumption 
is equivalent to the restriction to small Mach numbers and small 
Helmholtz numbers, where the Helmholtz number is defined as the 
ratio of the typical length of the problem to the wavelength of 
the sound waves. For our problem, the latter condition means, 
that the interaction region close to the lip should be small 
compared to the acoustic wavelength. It will turn out, that this 
region has a dimension of the order U Q /f (Vfhere U Q is the mean 
flow velocity and f the sound frequency) . Consequently, we should 
have U 0 /f« a Q /f, where a Q /f is the acoustic wavelength. This is 


equivalent to M « U 0 /n 0 * 1 < l.e., again the condition of sund t 
Mach number. The seventh assumption. the restriction to an infi' 
nitesimally thin shear layer will limit the validity range of 
the theory to the case where the shear layer thickness is small 
compared to the wavelength of the instability waves. In other 
words, the Strouhal number fO/u o should be small. O is the 
momentum thickness of the shear layer. 


After all these simplifications we end up with the linearized 
Euler equation and the continuity equation in the following form: 


( 1 ) 


(2) 

(3) 


-iusU + U 
-i<i>v + D 


an 

ax 

3v 

ax 


+ v 


DU 


3y 


+ 


1 Dp 
J 3X 


+ 12 2 

o ay 


o 


o 


D U 

Dx 



o 


where U is the mean flow velocity and u, v, p are fluctuating 
quantities, all proportional to e"^ wt . 

'fhe classical approach would be to fulfill the boundary condi- 
tions at both sides of the shear layer, which means that both the 
pressure p and the displacement h should be equal there. We 
focus first on the pressure condition. Instead of taking p we 
may an well choose Dp/Dx. Above the shear layer (index "1") we 
have U * o. Equation (1)then reads 

1 ap 1 

(4) - iuUl +-3^- - o, 

U Q and 


Below the shear layer (index "2") we have U » 

- au 2 1 0 P2 

(5) -i»u 2 + + = o. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


\ 

•f 


•a 



If we subtract: eq. (5) from eq. (4) we obtain 


( 6 ) 


u 2 + i 


u 3x 


_i_ 9 (Pi-p 2 ) 

pit) 3x 
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For x > o we should have P^~P 2 B Ap, 2 = °* 

This is a first equation connecting the u-velocities at both 
sides of the shear layer. 


Now we consider the condition of equal displacements. The dis- 
placement h and the velocity v are connected in the following 
way 

(7) v - + 0 |5 • 

In the region above the shear layer we have with 0 = o 

(8) v 1 «'■ -iuth; 

and below the shear layer we obtain 

(9) V 2 - -l»h t B 0 15 . 

We can insert h from eq. (8) and find 
5 o 3v l 

do) v 2 - Vl + i — ' 

which is the second connecting equation between the velocities 
at both sides of the shear layer at y = io, However, equations 
(6) and (10) are valid for two different velocity components. If 
we want to find solutions for v and u we have to generate (i) two 
equations for the two unknowns v^ and v 2 as well as (ii) again 
two equations for the two unknowns u^ and u 2 » 

(i) A second equation for v (| and v 2 can be produced in two dif- 
ferent ways. In [12] it is shown how eq. (6) can be con- 
verted into an equation for v using a source distribution 
approach. We will follow, however, another procedure (which 
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iB also given in tl 2 .1 ) which will provide more insight into 
the physics; this latter approach we will call the pressure 
gradient approach . 

(ii) To create a second equation for u 1 and u 2 we have to convert 
eq. (10) from an equation for v into an equation for u. This 
will be done using a vortlcity distribution approach . 

2. 1 . Differential equation for v 


The aim of this section is to create a second equation for v^ and 
v 2 out of a consideration of the pressure field and its gradients 
at the shear layer. Then we will derive a differential equation 
for v 1 only and solve it. 


We start out by taking the x-derivative of the first Euler 
equation (1). Subsequently, we take the y-derivative of the 
second Euler equation (2) . Both derivatives are added and some 
terms are eliminated using the continuity equation (3) . We end 
up with 


(ID 


V 2 p 


-20 


8 V _ 3_U 

3 x 3y 


In our model, the mean velocity profile jumps from 0 = o for 
positive y to U = U o for negative y. This can be represented in 
terms of a Heaviside step function H 

(12) 5 = U Q • H ( _ y) . 

The derivative of a step function is a delta function and we 
obtain from eq. (11) and (12) 

(13) V 2 p = 2pU • • 6 , , . 

■ o 3x (y) 

Eq. (13) is a nonhomogeneous Laplacean equation with a line 
source of varying strength in the shear layer. It should be 
stressed here, that the sources in the shear layer are pressure 
sources and not sources of fluid? the latter would violate the 
continuity equation (3). Since there might be some confusion 


15 


ORIGINAL PAGE IS 
OF POOR QUALITY 


about the structure of these pressure sources we will give, as 
an example, the detailed pressure source structure of an ampli- 
fied instability wave in appendix B. It will be shown there that 
eg. (13) can be integrated to produce the whole pressure field. 
For the further progress of our calculation, however, a detailed 
knowledge of the pressure sources is not necessary. 

At the surface of the semi-infinite plate we have v and 3v/Dx 
equal to zero. Consequently, the pressure source strength is 
zero on the plate surface (see eq. (13)). The only other location 
where V 3 p is nonzero is at the location of the exterior pulsat- 
ing source (see Fig. 2). The basic idea of this approach is, 
that the pressure distribution in the whole field can be split 
into two contributions: 

(i) a pressure field which is symmetric with respect to the 
shear layer and which is created by the pressure source 
distribution in the shear layer itself. 

(ii) a pressure field which is created by the exterior forcing, 
e.g., a pulsating source. The pressure fluctuations of this 
contribution are transmitted through the shear layer. The 
pressure gradient of this contribution is continuous through 
the shear layer and therefore it is antisymmetric close to 
the shear layer. 

As a result of this splitting process we have 


Pi = P lB + P 1f ; P 2 = P 2s + *2£ 1 

<14)<j v, - v 1s + v 1£ ; v 2 = v 2s ♦ v 2f 

U 1 = u 1s + U 1fi u 2 = U 2s + u 2f 


The index s stands for shear layer and the index f labels the 
exterior forcing. The boundary conditions at both sides of the 
shear layer have to be fulfilled by the summations of the 
individual constituents, e.g., by v 1 and v 2 ; as before. On the 
other hand, we have some new informations: Since the induced 
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field of the pressure sources in the shear layer is symmetrical 
(it is created by Sources of symmetrical directivity in a field 
with symmetrical boundary conditions) we obtain: 


(15) 


a Pls 

oy 


a P2s 

3y 


at 


y « ±o 


and for the continuous pressure of the exteriour forcing 


(16) 


3p 


If 


oy 


3 P 2 f 
‘ ay 


at y 


±o 


These conditions for the pressure gradients in y-direction can 
be inserted into the 2nd Euler equation (2), which gives 


(17) 

v 2s + 1 

^o 

U) 

3v 2s 

3X 

“ “ v 1s 

(18) 

V 2f + 1 

w 

3v 2f 

3x 

s v 1f 


Equations (17) and (18) can be added, using eq. (14) 

0 3v, 

(19) v„ + i + v. - 2v 


If ’ 


This is the desired second equation for v 1 and v 2 - The velocity 
v.j£ is not an unknown quantity, it is the Velocity which is 
created by the exterior forcing without the mean flow being 
present, but in the presence of the semi-infinite plate. 
Equation (10) can be inserted into eq. (19) to obtain a non- 
homogenous differential equation 


( 20 ) 


2v. 


+ 2i 


£a fn . f!ii 

w dx t id J dx 2 



The homogeneous solutions of this equation are found easily 
(21) cp 1 = C^l*, cp 2 = C 2 e A 2 x 
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C 1 and C 2 are arbitrary constants with the dimension of a velo- 
city. and Cp 2 ar ® the well-known spatial stability waves for 
an infinitesimally thin shear layer, extended from -» to +», 
which we will call the Helmholtz solutions . The complete solution 
can be obtained in a straightforward way 


e ^1 X v^^dx + 


e ^2 x v^dx 


This is the general solution for v^ . Once v^ has been found, v 2 
can be determined with eq, (10). The lower boundaries and 
G 2 will produce terms of the Helmholtz solution type. The actual 
value of G^ and G 2 depends on the boundary conditions! this will 
be discussed in detail for our configuration (shown in Fig. 2). 
However, before we are going to do this, we have to consider the 
structure of v^£, the fluctuating vertical velocity at x > o 
and y = o which is created by the exterior excitation. 


The general solution (23) for v, was obtained by splitting the 
pressure field into two constituents p^, the forcing field and 
P s the field induced by the shear layer Wave motion. The field 
p^ is continuous through the shear layer. It can be seen, e.g., 
from equation (2), that the velocities cannot be continuous 
through the shear layer, in other words: v^£ ^ v 2 £. Therefore, 
the appropriate approach is to calculate the pressure field 
first , which follows from V p = o outside the source location and 
the shear layer. After splitting it into the two constituents we 
have (i) the excitation pressure field p£ with a singularity at 
the location of the exciting source (x q , y Q ) (see Fig. 2) and 
(ii) the shear layer pressure field p s with a pressure source 
distribution in the shear layer (x > o, y = q). As we know from 
potential theory, analytic functions of z = x ± iy are solutions 
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of tho haplacean equation foe the pressure. Therefore, wo can use 
all known tools of conformal mapping to calculate* o.g«* tho 
pressure field p f . in tho fluid at rent (0 * o in region "I" of 
Fig* 2) there in a simple relation between the velocities and 
the pressure. This can ho obtained from the two '.Sul or equations 
(I) and (2), 


(24) 


U. « 


1 ilp 
r u’ i x 


i v. 


. .1 ip 
cw ay 


We recognise* that, in a fluid at rest, tho pressure is pro- 
portional to the potential function which is defined as 

( 2 r ») u a At ; v o jt * 

This allows us to use conventional potential theory to calculate* 
say, v | p for a source or a dipole in conjunction with a semi- 
infinite plate. Tills has boon done in [1*2] and we will give utily 
the results hero. For a pulsating two-dimensional source, as 
shown in Fig. 2 wo have 


(2G) 


'If 

(y°o) 


> x 


0 ?o * * " jj*o ” x n^ 

• Tit ' Tx - ) V 1 y 3 “ “ 


with 


(27) 


r 'j | x ^ \ f • 

O O 3 0 


x Q and y 0 are the horizontal and vortical distances* roapoetivo- 
ly, of tho source from tho plate edge and 0 is the source 
strength, say, in m J /s, Clearly, Q is fluctuating as 


(20) 0 * Q 0 0 


-iwt 


For large distances r » x we have simply 


(29) v. „ » - -i 


''lf---= , £V§- „ - V r„). 

u 
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in 112) the induced field of a dipole is also given, we will not 
give hero the full equations of the general ease of an arbitrary 
location of the dipole. There is, however, a simple case which is 
worth mentioning. Assume that the dipole is far away from the 
edge, i,o., r >> x, and located above the end of the plate so 
that y 0 » x Q (see FifbJ) . 


-X. \ 




\ 


dipole 


Vo 


i 


A 


plate shear layer 


Pig, 3 Dipole ns excitation of the shear 
layer. 


The direction of the dipole is \V . In this case, wo have, 
according to t!2): 


(30) 



* eos ($- ) , 

■' ,Iy o 


where D is the dipole strength. As we see, 


at » 4t 


v ( £ becomos a maximum 


The direction of maximum efficiency of the 


dipole clearly changes with its position relative to the semi- 
infinite plate (see also [12]). For a position close to the shear 
layer, wo will have, as expectod, \^ mnx " 90° . If wo use a vibrat- 
ing cylinder to create a dipolo, we have for the dipole strength 
O a u ■, SiitR* where R is the cylinder radius and u„ is the velocity 

L* ,*4.1, L, 

u c au o‘° ot ’ ^ho Similar expressions can bo given 
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for a vibrating ribbon in the fluid at rest tl2l.*) The general 
finding is that in almost all conceivable cases we create a field 
close to the plate edge where v 1 f « 1//x. There are only a few 
singular cases where we create a stagnation point at the plate 
edge, like the "zero" position of the dipole in eq. (30) , In 
section 4 it will be shown, that this region close to the edge 
where v^£ * 1/*dT is in fact the relevant one for the whole inter- 
action process. 

As we will show now, v^£ “ 1/^x corresponds to a parabolic 

pressure field around the edge. This fact is of particular signi- 
ficance for the choice of a suitable reference quantity for the 
excitation in real experimental situations. Since it is often im- 
possible to track down the whole pressure field in a practical 
situation, we are looking for measurable quantities relevant to 
the interaction process. 

To model a parabolic pressure field around the semi-infinite 
plate, we have to keep in mind, that analytic functions of z=x±iy 
are solutions of V*p f = o and that we can use conformal mapping 
to fulfill the boundary conditions at the semi-infinite plate. We 
start with the consideration of a plane pressure field (see Fig. 4 ) 


Z n - plane Z h -plane 



Fig. 4 Plane pressure field, parabolic pressure field, 

*) If the vibrating ribbon operates within the flow regime at 
y < o, the generated pressure field is, however, completely dif- 
ferent, due to the interaction with the flow. 
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This pressure field in the z a -plane <z fl * x a + iy fl ) has vertical 
lines as lines of constant pressure (we suppress the time depen- 
dence e"‘ 1 ' wt for convenience) . This would be the pressure field 
of a horizontal velocity fluctuation u constant throughout the 
plane. In the z a ~plane we have p *» A • X a * If we apply the con- 
formal mapping function = z fl 2 we have for a point 
z a = x Q + iy a = r^ 1 ” a conversion to the point z fa = x b +iy b “ 

r.e"^ w z 2 s- r 2 e 2ia (see Fig. 4). 
baa 

For the pressure we obtain: 

(31) p = A • x g = A • Re(z a ) = A • Re(*^Zj^) <* A ’ • cos , 

For 8 = o (above the semi-infinite plate) we have p 1 = Av£^" an<J 
for 8 = 2 n we have p 2 => -A/x^* (below the semi-infinite plate). 

At the location 1 of the reference microphones we have 

(32) p.,-p 2 “ A P 12 = 2A , 



With this equation we can eliminate A and find 
Ap 12 / — 8 


(33) 


P = 


2fr 


cos ^ . 


To obtain v 1<: we take the derivative at y = o and x < o, which 
i 1 i } D ay 

is equivalent to - — at 8 - n. 


(34) 


13] E = . Ap 12 

r 38 ( 8 = it ) 4 


Since our plate (see Fig. 2) should be extended from x = -“ to 
x = o, and since we want to have our result in Cartesian coordi- 
nates, we write 


AP 


12 


(35) Ty (x >o, y = o) = + 

Using eq. (24) we can relate this expression to v^ f 


(36) 


If 


1 


^12 
4 pw/l 
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This is the desired relation between a quantity easily mensurable 
such ns fcne pressure difference A p. 2 at: the plate taken at a 
distance 1 from the lip, This quantity, &p <2 , CG htains ho con- 
tributions from the shear layer. We feel, therefore, that A P| 2 
is the most suitable reference quantity for the excitation of 
the shear layer. 


Since we have with oq, (36) a distribution of the excitation 
field V | £ Which is valid in almost all conceivable cases, wo can 
proceed to solve eq. (23), the differential equation for the 
shear layer motion. After inserting v.. f (eq. (36)) into eq. (23) 
we have 


(37) v 1 * 



•M* 0 >T 


The lower boundaries s. and G a will produce contributions to the 
solution of the Ho Imho ifest- type, Since v, must be aero for x < o 
those are truncated Uelmholtss-solutiona starting at x «a q, fn 
appendix C wo prove that such solutions cannot fulfill the con- 
dition of equal pressure at both aides of the shear layer and 
that there As no combination of G, and g 2 fulfilling this con- 
dition, ’I'hia calculation in the appendix is based on calculating 
vs from v, using the condition of equal deformation of the shear 
layers (eq. (10)), Subsequently the corresponding u., and u 2 
distributions are calculated using a source distribution approach, 
Finally, using the first Hu lor equation (eg, ())), ,1 p ( /? x and 
ap^/Ax are compared. Since ;>p.,/;\x and ap^/ax are incompatible 
under all conceivable situations, we can conclude, that there 
are no unforced (free) vibrations of the shear layer , This is, 
by the way, one of the essential results of the present investi- 
gation, 

Since wo have G 1 » 6j * o, we are loft with a forced solution 
only, for which we have to prove that it fulfills the conditions 
of equal deformation and equal pressure at both sides of the 
shear layer. This is done in several stops* 
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(i) v 1 is calculated in an unambiguous way from oq. (37) with 
0^ * 0 2 * o. subsequently, v 2 con bo determined in an un- 
ambiguous way using oq. (IQ) , 

(ii) by a separate procedure (section 2.2.) analytic solutions 
for u 1 and u 2 are generated. We end up with a differential 
equation for very similar to eq. (37) , The free constants 
there are determined by the compatibility with the and 

v 2 solutions in the far field. 


(iii) since we found no direct approach to relate the analytical- 
ly obtained to and u 2 to v,, in the near field, this 

was done using a numerical procedure, it is a by-product 
of the calculation of u 2 in the whole field (not only for 
y * o) which had to be carried out for a comparison between 
theory and experiment. 


(Iv) it is shown, that the analytic solution for u^ and u„ loads 
to zero pressure difference at both sides of the shear 
layer ( x > 0 ) and also to the expected pressure difference 
of the excitation field at the semi-infinite plate (x « o) . 


As we sea, all ambiguities and, hopefully, all doubts are re- 
moved by this chain of procedures. The first step (i) is then 
to evaluate eq. (37) with * G 2 * o. After some suitable sub- 
stitutions (see (12) for details) and introducing the error 

function defined as 
z 


(38) erf a 

» J. 

■ 

rir 

wo have 

0 

(39) v^ * 

AP 12 . 

PtJ‘ r l 

O 

and with eq 

. (10) 

( 40) v 2 * 

ip , 2 



i »"tf 
1 


X |X XoX 

e 1 A""r: _ Q 1 




erf A ^ x 


erf a 2 x 


A 

A 


' \ 1 X 

— — — erf '‘TTx + - erf ^x 

./c 1 a ; 2 

1 •* 
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( 22 ) 


V B (i + 1)) ^rj ® (i " 1 ) , 

' tl * r? 


It is useful to define dimensionless quantities as follows: 
A dimensionless wave number 

(■II) ^ “ i + 1 1 ■ i - I , 

a dimonaionloas distance 


(. 12 ) 


x “ 


■dK 


which will bo used in the same way for y; and a dimonsionlosa 
velocity *> 


(43) 


v 


v 

a e: ** 


• »' • aj, 

Ap 12 


which depends on quantities which can be measured easily. (The 
same form will be also applied to u. After introduction of those 
dimensionless quantities we have* 




(44) 


V, a 

I 4 


i * n | Q V 


Lb 


5?*" erf y V. x 


( 4 $ ) | 53 ** 


' \st 
e_ 1 


~=S“ erf/L,x 
$2* 


«v j^f***^» ^ ‘V 2 « rr™®* 

Q££ Lit Hh ***^^ i *' A 61 *f V 


* V> 


L. 




Pin . 5 shows a plot of eqs. (44) and (45) where and v^, have 
been split into their real and imaginary parts, respectively . 


*) Based on dimensional considerations, similar dimensionless 

quantities had boon defined in M2]j the relation to these is: 
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Tables of numerical values of v 1 and v 2 will be shown in appendix 
D« Subsequently, we will investigate the behaviour for small and 
large because this provides some additional insight into the 
interaction process. 


For small x wo take advantage of the expansion given in f 3 1 J 


( 46 ) 


erf z * 


E 2 n 

n*o l" • 3 . . . (2n+l 


,2n + 1 


wc obtain for small x 

(47) ^ « 5 1,5 i v 2 4T’ 

which shows, that thore is no singularity of v., or v 2 anymore 
near the trailing edge, in addition, it is shown by eq. (47) , 
that the flow loaves the trailing edge without a jump in the 
slope. This is equivalent to the presence of a Kutta condition . 


A question of particular interest is the far field behaviour of 
equations (44) and (45), By virtue of the far field expansion of 
tho function erf z , given in (31) 

2 » 

(48) /n za 7 orfc 2 l 1 + E (~1) m - — , 

2 •» » m * 1 ( 2 z * ) 


with erfc z* 1 - erf z we obtain for large 2 , where the 
summation term wan be neglected: 


(49) v 1 

(x + * 


1 Sir <3 *1 * __i__ 

4 ‘ " 4/jf 


(50) V, 

4 (x + <° ) 



•t; 


i 



with X,. « i + 1 . It is evident, that the field has separated into 
an amplified Helmholtz-type wave (1st term of eqs. (49) and (50), 
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respectively) and the contribution of the exterior excitation 

which is the second terra in both equations. 'Phis latter con- 
tribution is transmitted through the shear layer without any 
change. If we call the first term in the equations v f } (where H 
stands for Helmholtz solution), we may write equations (49) and 
(50) 


(51) 


(X^“) 


V 1H + V 1f 


' (X+“) 


V 2H + V 1f 


It is interesting to see where this far field solution emerges 
i.e., beyond which distance x amplified instability wave and 
excitation can be split. In order to investigate tliis, we sub- 
tract the asymptotic instability wave from the complete solution 

? ? 

< 5 2 > v 1D " v 1 - v 1H * v 1f ; v 2D = v 2 - v 2H - v 2f . 

The differences v 1D and v 2D are plotted in Fig. 6 . 

It is evident from Fig. 6 that the far field solution is valid 
beyond x * 4 t 6, that is beyond about 2/3*1 wavelength of the 
instability waves. This means, that the basic interaction takes 
place within ^ wavelength of the instability waves. This latter 
finding, that the wave field can be separated Into a predictable 
instability wave and the exterior excitation field (sound field) 
is one of the essential findings of this paper and it supports 
similar attempts made previously on a more heuristic basis (32] . 
This splitting is therefore admissible at distances greater than 
about 1/2*1 wavelengths of the Instability waves. Finally, we 
provide also a simple formula for the magnitude (modulus) of the 
excited instability wave: 

(53) I v., 1 
H 


= Iv, 


/7 x 

r ■ ■ rr-nriT * 0 

^ ^ for x - 4 * 6 
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(54) 


r. «ri^m:^:»r3s=ai!ncsaic jm,-;nr*r*ir» t. ;*fc j*s »« 


lv 1 I * I Vg I « 

'll ^ii n* f u 


M 

0 % 


4-'V2 


Numerical valuos of v. and v, nro given in appendix D, 


Since v-velocity components are difficult to measure, we will, 
in the subsequent sections, both calculate u^ and u^ at y « 1 o 
(analytic solutions) and will provide also numerically computed 
values of u 2 outside the shear layer in tha potential field. 

2.2. Differential equation for u 


The procedure to obtain a differential equation for u is slightly 
less elegant than the one for v. Wo restrict ourselves, in 
addition, to the "parabolic" excitation, for which wq obtained 
the special solutions, oqs. (4 4) and (45) in the preceding 
section. As wo have shown before (see also section 4) the "para- 
bolic" excitation ease is valid in almost all conceivable 
situations. In addition, when dealing with u-velocitios, we do 
not worry about causality? and if wo encounter ambiguities, wo 
can remove those by linking the u-solutions to the unambiguously 
known v-solutions. If wo create, e.g., a parabolic pressure 
distribution around the semi-infinite plato it is irrelevant, 
therefore, whether this has been created by a singularity distri- 
bution at y « 'o or by a source at a large distance. As a con- 
sequence thereof, we do not split the field here into two con- 
tributions from the exterior source and the shear layer itself, 
respectively. 

The first task in this section is to construct an equation for a 
out of the equation for equal displacement 

no) v, - v, + iSa U- . 

In (12 J and in appendix C it Is shown how n v-velocity condition 
can be croatod out of a u-velocity condition using a source 
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distribution approach* Hero wo will do the opposite, i.o., wo 
will go nor ft to a u-veloeity condition out of ft v-velocity con- 
d it ion using a vorticity distribution approach. 



ip 2 

Fig. 7 Induction of a vorticity 

distribution on a straight lino 


Wo. consider the velocities which are induced by a vorticity 
distribution on a straight lino between b., and b 2 (see Fig . 7 ) . 
Tho induced velocities u and v at; a point P^ above a vorticity 
distribution -rtf.) on the x-axis aro: 


(55) < 


u a - JL V(i) ’ Yl 

2 « J (x-f .) 5 'ty, 3 ' dr *' 
b l 


2 it 


Ik] • (x-f.) 


dt. 


if the point P. ( approaches the x-axls (+y j *■ o) , the u-Gomponent 
becomes u « ~y(5,)/2. This fact can be verified if one analyzes 
the distribution yU) into a series of small stops of locally 
constant -y U) end evaluates eg. (55). One finds that in the case 
of a vorticity distribution on a straight line no induced 
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u-component contributions can come from vorticity elements out- 
side 6 (at +y 1 **• o). By inspection of Fig. 7 this is also obvious. 

For a point P 2 approaching the x-axis from below (-y 2 •*■ o) we 
have u = +yU)/ 2. With these considerations we obtain: 


(56) < 


(y=+o) 


(y=-o) 


b2 

1 f u ( e ) 


ir J X-5 
b 1 


4 

IT 


u(F,) 

X~E 


dE 


at 


Inserting eq. (56) into eq. (10) and adjusting the limits b 1 and 
b 2 yields 


+« 


(57) 


' u 2 
x 


U-U) 


d? = - 


u 1 U) 


d? - i 


U -S U 1 

O ■! 1 

U) Sx j X' 


u. (e) 


■5 


d£ . 


Since we expect exponentially growing instability waves in the 
downstream direction we will run into difficulties with the 
upper limit of the integrals in eq. (57). To avoid this we will 
use a simple trick: we will subtract from the forced instability 
solutions the amplified Helmholtz solution u H which has asympto- 
tically the same magnitude as the forced instability wave. 
Consequently, we write 


+« 


+ W 


(58) 


(u 0 -u 2H ) 


• iu 2 - 
x- 


d5 + 


<u r u iH> 

x-e 


u. 


d£ + i 


u) 3x 


<u r u iH> 

x-c 


d? = o, 


The evaluation of the third term in this equation uses 
integration by parts 


(59) 


_3_ 

8X 


f <U 1 ( E ) ~ U 1H ( E) 
X-E 


<U 1 ( E ) ~ U 1 H ( 5 ) * 

(x-U s 
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« W T* 


The fir tit (integrated) form on the right hand aide becomes aero 
if (i) x remains finite, (ii) the difference (u,-u )n ) does not 
grow faster than t, for large negative and positive values of f* 
That the latter condition is fulfilled will ho evident from the 
solution we are going to obtain*, so wo have 


(f>0) 


-i <n 


(U 2 -U 




211 ' 


x-f. 


clt * 


x-r. 


cH + i 


;o 

r u> 


f 


r&v 


,u tn> 


X-t 


at , 


The integrals in uq. (60) each posses the same denominator and 
the same boundaries. No have a condition at least aa restrictive 
if we write 


( 61 ) 


VI r. "* U ■ 1 


* u. 


u. 


in 


i 


jo 


a 

'US 


(U, 


* u ui> 


« 0 


whore wo have replaced h by x. Since wo know that eg. (61) is 
fulfilled by the UolmholU solution alone wo can write 


;-nrr».*:i*s -m -Mu^taotSisSX' "~:rr xv*#*#®*? r"*3«*atvir-3msi 



it is valid for all x since eq. (10) was valid for all x. This 
is the desired second condition for u. ( and u,* if we insert this 
equation into eq, (6) we obtain a differential equation for u^ 
alone: 



No can obtain the pressure at both sides of a semi- infinite 
plate (located at x < o) from eq, (33), after adjustment of 
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J 


(y«-o) 2 rl 

This can be used to find the bight hand side of esq. (G3) 


l) du. u 2 d u 1 ,, i Ap .. 0 
(65) 2tt,+2i -7^“ (-2) — -yl » -4™ • ~4 

1 u ' dx * ax 2 /=* 2 PW /r 


for x > o the right hand side of oq. (65) becomes zero, because, 
by definition, the pressure difference at both aides of the 
shear layer should bo zero* In complete analogy to the procedure 
to solve the differential equation for v in the preceding section 
wo obtain: 


( 66 ) 


*^12 

4i'U 0 /i 


X X 


0^ f 

°"" s , 

lx - o’ V 2 x [ 

tf *2* 1 

dx 

L J 


J 

• -"x “•* 


We should keep in mind, as mentioned before, that thero exists 
a contribution to the integrals only for x <o, Consequently, we 
have for x > o only contributions of the Helmholtz instability 
wave typo , created by the lower boundaries of the integrals 
and G f j in oq. (66). 

?n unambiguous way to find the magnitude of these Helmholtz-type 

contributions is to rolato them to the asymptotic magnitude of 

the v instability waves. Those v instability waves also behave 

like Helmholtz- type solutions for distances ££ £ 4 * 6. We 

llo 

should find out, therefore, what the coupling between the u and 
the v-Gomponents in the Helmholtz instability waves is. The 
easiest way to do this is to start ovjt with the potential of 
the induced velocity field of these waves? 
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we take the derivatives 


v 1 * iA^ 2 Ao x 1 ,2 i x+i y> 




v 2 * -i.\ 1>a Be X 1/2 <x * ;L y> 




Therefore we have u^ * -iv^ and u 2 * iv 2 at y ■ o. With the 
solutions for v^ (eq» (44)) and v 2 (eg* (45))/ we obtain (with 
u nondimensiQnalized as in eg. (43) )t 


- A r“e A I 5 e' V 2 x 1 
u i * T — “ ~-J 
’ *1 ’A 2 j 


r A .*\X 

. r tt e 1 ^ e 2 

4 /*—* /zr* 

^ / *2 J 


for x > o 


u • p /idU 1 


O . , MX ? T 

, * * — f -\ * i + i; « i - 1, 

u A 

u o 


Analytic solutions of this type had not been published before/ 
but they were obtained already with different mathematical 
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methods by W. Mfihring 133] and M.S. Howe 1 3 4 3 . To find solutions 
for x < o we have to proceed analogously to section 2.1, The 
only difference is, that we now have additional Helmholtz-type 
solutions in contrast to the situation for the component v. After 
some algebra we end up with 


U 1 n T 


Pe A 1 x „ „ e X 2 x 

erfc/^* - — 
u ai *^2 


ij r ercc ■v 

L ~ y\ * 


cfc ^ 2 X J 
erfc 


x has to be inserted into these equations with it-s negative sign. 
The abbreviations are identical to those used in eqs (69} , (44) 
and (45). Since erfc is defined as erfc z = 1-erf z, it is ob- 
vious that equations (69) and (70) merge into each other at 
x = o. If one expands eq. (70) for large negative x, one can 
show that the wave structure of the instability Waves vanishes 
and the prevailing terms are proportional to u « 1//-X. In 
addition, it can be shown using the first Euler equation (1), 
that eq. (70) produces the correct pressure difference at both 
sides of the semi-infinite plate. 

For large positive values we have from eq. (69) 


lu I - iu [ - APl2 

- ^ c 

tax 


j 

1 ' 1 2 

4.V2 

Dp 


.. _ i 

$ 

for x Z U o /f 




i 

the same result as eq. 

(54) for 

v at 

large x 

1 

t 

, The only 


difference is, however, that the excitation velocity has no com- 
ponent in the u-direction for x > o. 

There is yet another less obvious feature of the exact solu- 
tion, eq. (71), near x = y = o. If one wants to verify eq. (71) 
experimentally, one might take, e.g., data of the velocity u 2 






in the region outside the shear layer and would extrapolate these 
data into the shear layer. However, this would lead to erroneous 
results near x » y « o (plate edge). The reason is, that the 
gradient of a curve u 2 (y) for constant x would be infinite at 
x * y • o, which hardly allows any sensible extrapolation 
towards y o from data at y *j> o. The reason for this is the 
following. In the potential flow region outside the shear layer 
the vorticity should be zero, i.e. , 3u 2 /3y - 3v 2 /3x = o. Since 
v 2 “ /x for small X, Du 2 /3y must become infinite (but, certain- 
ly, not u 2 itself), Thus, the exact solution for u 2 (eq. (69)) 
is only of limited value if it comes to comparing it with experi- 
mental data. This is the reason why the u 2 velocity field for 
finite y has to be calculated numerically. After this will have 
been done in the next section, and the Induced u ? -fleld is 
available, the discussion in this paragraph will also become 
more convincing. 

3. Field calculations 


For a comparison with experiments it is necessary to calculate 
velocities also outside the shear layer. The quantity which can 
be measured most easily and accurately is the velocity u,, in 
the flow region y < o. A source distribution approach is used 
to compute u 2 from v 2 , which is given for all x. Actual data 
of V 2 , i.e., its real and imaginary part can be found in 
Appendix D and curves can be seen in Figs. 5 and 6 in section 
2.1, The source distribution approach we use here is identical 
to that used in [12] or in Appendix C; there is also a formal 
similarity to the vorticity approach used in section 2.2. 

We consider first a source distribution q (£) on the x-axis 
(see Fig, 8 ) 
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Pig. 8 induction of a source distribution. 


The induced u- velocity component at a point p 1 above a lino 
source distribution q(f) i3 


(72) 


u * 


92 

_L f q(f.) * (x - r.) 

2w > (x, - u a + y, 2 
g, 1 1 



Pox* a point Pg below the line source distribution we end up with 
the same result. If the point P 1 approaches the x-axis (+y^ +o) 
the v- component becomes +q {«;)/?.. For the case of a source distri 
bution on a straight lino no induced v-component contributions 
can come from sources outside f at y_| ►o, So q (C) can be re- 
placed by 2v for y > o. For negative y we obtain negative induced 
v-components from positive source strengths q (£) . Thus we havo 
for negative y 


(73) u 2 - 



CJ 1 


v.j U) • (x,,-E) 

- A — 2 — i — - • 

(x 2 -o +y 2 


This is the type of equation which provides the relation between 


v 2 and u 2 , It will ba evaluated numerically, Formally, we would 
h&ve to integrate from g ( * -« to g 2 # + *, However * we again ©n« 
gaunter difficulties with the upper limit* because growa 
exponentially. This difficulty will be circumvented in the same 
way as wo did it bo Core. We sub tract; the exponentially growing 
asymptote v 2R to obtain an integrabl e function. The difference 
V 2D * v 2 ~ v 2ji iB on ly proportional to l/v'x for large positive 
x and becomes aero exponentially at largo negative X. After 
having carried out the integration we add the induced field off 
the asymptote v 2 , } (see eq. (68)) to obtain the complete u 2 -dia- 
tribution. This whole procedure lias to be carried out separately 
for the real and Imaginary part of v 2 to obtain the real and 
imaginary part of tu* 

The integration of u 2 is then split into three regimes*. 

(i) For largo values of x the real part v 2DR of v 2Q becomes 
aero and the imaginary part v 2oI of \' 2U is proportional 
to t//x (see Fig. G). Thus , for values between, say, 
x « G and x * w© will use the asymptotic form off v«, D j, 
insert it into oq. (73) .and will carry out the integration 
in closed form. This will take care off the contribution to 
the integral (73) for large values of x (or C) . 

<U) For small positive x the real part v 2r of v 2 has a para- 
bolic behaviour with v 2R « M, Caused by this* the slope 
of v 2 p R is infinite for X * +o. Since this can lend to 
numerical problems wo will also integrate the region be- 
tween k « o and x * o.1 in closed form, 

(ill) The remainder of the integral from, say, x <* -6 (whore 
v 2 * o) to X ■ o and from x * o» 1 to x * +6, where the 
asymptotic solution takes over, is calculated using a step 
integration procedure with intervals of Ax * o.T. The stop 
Integration elements are trapezoidal pieces where the 
upper and the lower limits have the same values as the 
function v 2dr (or v 2DJ ) . This provides a much better 
accuracy than ordinary rectangular steps would do. 
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For each of these Integration elements (i) * (lit) we have to 
produce analytical solutions, We start with the far field in- 
fluence (i). For large values of x we have for v 2DI from eg. (50) i 


(74) 


° 2DI * " 7r* 

(far field) 


For u 2 , x and e we use the nondimensionalized form to obtain 


(75) 


G 2di ** + T* 
(far field) 


CM 


J 

h 


(x-l) 

<x-t) 2 -y 2 




We will not find this type of integral in the usual integral 
tables. However, it can be simplified by complex fraction de- 
composition Using the following relation 


(76) 


JhL 


*v s >v 


1 

2 




with 

(77) z « jt + iy » z = x - iy , 
Thus we have for (75) 
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(78) 


2DI 


1 

8* 


d5 


(far field) ^ 


/| (z-l) 


F? J 


de 


(5-D 


Integrals of this kind can be found, e.g., in (38). The total 
integral turns out, as it should, to be completely real. We 
obtain after some algebra: 
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il5 


(79) u 


staffs 

, r 1 ,» r+c 1 - 2 / 5 Tr cos ^/2 

„ —L- i • cos t • In <- ^ -— 

201 4tj / r »— 2 2 r+5 1 + 2/^r cos\>/2 

A 2/f7t sin £/2 a- 

+ sin ^ • arctan < ) -it sin j 


) + 


with r <* /ST"* y 2 and $ * arctan (y/x) , and for 
in the order of 6. 

The second <ii) particular element of the integration is tackled 
in a similar way. The expansion for v 2R (the real part of v 2 ) 
yields for small x 


(80) 


v 2r - -/£, 


Therefore, we have to integrate 
? 2 


(81) 


U 2R * + J 
(small x) 


<2-^ as 
is-Si 2 *? 2 


for f, 2 in the order of o. 1 . we proceed in a way very similar to 
that shown above and end up after some intermediate calculations 


(82) 


, /- r i A r+ h~ 2/ h* COB 

“2R ’ f [- 5 ' “» I ’ ^ — ) * 


(small x) 


r+l„+2/fTr cos $/2 


4 


+ sin 2 ’ arctan ( 


2 /FTr sin $/2 


2/r, 


r-£. 


) 


f ~i 4 r 
• u 


with r and **» defined as before. 


This result has most terms in common with eg. (79) which 
suggests, that the same computer subroutines can be used. 


<11 






The third (iii) element is the trapezoidal element mentioned 
before (see Pig, 9) 


/ V2DR \ 

(or V2D1) 
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Pig. 9 Trapezoidai velocity element as inducing source element* 
v * * \ • source strength • 

The local distribution of Vj D (in the computation specified as 
V 2DR a:id V 2DI> is 


(83) ^2D = °a + m {r > - c a* 


(84 ) m 


v b " v a s v b ~ v a 


The integral 


(85) AUr 


(v a +m(5-5 a )) • (x-fj 


1r (x-fj 2 +y 2 


can be solved analytically using the integral table in [35] . One 
Obtains after some intermediate calculations: 
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( 86 ) 



V\, 


l t»»X 

I y (am an ( -=-= ) -am an ( - ^ ) ) - 
1 V V 


V u 

-.tt— r-r 

V. 

b a 


a -x) ) 


In 


■n u 

1 b 





with 


W) 


(x 


,)“ * y* and i a “ 


(x - i.J“ (■ y“ v 


Thin oli' wont (HI) will bo uood In all all: tuitions whom v a|) la 
not acini or where the asymptotic solut Ions (I) and (ii) are not 

UlK’d. 


At tho ond of tho oompuf at tons, tho atiympiol lo solution 
Unslnq ihj. (08)) in added to obtain u s * u ao + u., n * ‘the details 
of tho computations are given in appendix li o# wall an tables 
with ooiiiputod values of u v , lloro , wo will show tho modulus of 
u., , because it. can bo measured In a at raiyhl forward way, in 
Fly. m a computed diagram of |u a ! as a function of y for 
yarlium x can bo soon, Wo have unod n logarithmic plot on tho 
vo at l on I axis of thin figure. Thin allows to 1 dent i fy rog ions 
of tho curves which show exponential behaviour? they will occur 
as straight linos, 

After inspection of Fin* 10 wo find that tho curve a for |u«j| at 
x > l show an exponent lal decay in the -y direction. With 
equations (68) and (71) wo can predict the induced velocities of 
the amplified Instability wave ; we obtain 

(88) liU * -'.il, . e x * e"l*l , 

•i • V a 

( X s 1 , Tj) 

This asymptotic behaviour is also given in tho dotted linos in 
the upper half of Fig. h, As we see, eq . (88) provides reason- 
able predictions already for x> 1.8. This Indicates clearly, that 


\ 

i 
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downstream of, say, x * 1*1.5 the instability waves dominate the 
induced u-field in the potential region outside the shear layer. 
It can be assumed, that this feature is exhibited more clearly 
in the u-component than in the v-component. This is due to the 
fact that the excitation has no u-component for positive x and 
small y (compare equation (69) for u and equations (49) and (50) 
for v) . 

One might argue, that there is also an asymptotic equation for 
fu,| at large negative x. In the region upstream of the plate 
edge, it is reasonable to assume that the excitation field do- 
minates. Since the excitation pressure field is known (see eq. 
(33), note that x is inverted there), the u-velocity field can 
be calculated using the first Euler equation. 

(1) -iwu + u + ! -|E to o. 

This is, with p known, a nonhomogeneous differential equation 
for u which can be solved in closed form and leads to a solution 
containing error functions of complex argument. This can be ex- 
panded for large x to provide a simple far field solution. We 
can produce, however, the same far field solution in a much 
simpler way. It is a reasonable assumption, that the gradients 
of the fluctuating u-flow become smaller with increasing distance 

from the plate edge. If, e.g., u is « l/Z^o^Su/ax will become 
3/2 

« 1/(-x ) which means that 3u/3x will decrease more rapidly 

with increasing (negative) x. Consequently, for large negative x 
the fluctuating flow will behave as if no mean flow were there. 
After neglecting the second term in eq. (1) we obtain a simple 
solution, i.e., the parabolic flow created in a fluid at rest 
by the pressure distribution, eq. (33). We end Up with 

(89) |u 2 | = ■ — i 1 - * cos arctan |?|). 

4 .Vjp+y" 7 x 

(x < -1) 

This asymptotic equation is also plotted (dotted lines) in the 
lower half of Fig. 10. It works quite well for x smaller than -1 . 
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We are left with a situation which can be considered to be ideal 
in computational fluid mechanics; The region where information 
relies on pure computation is limited and is imbedded into 
asymptotic solutions. How well the computation works is also 
shown by a comparison between numerically calculated values of 
u 2 at y = o and analytic values from eq. (69). 


In Table 1 these two sets of data are compared. The deviations 
are minimal in spite of the fact that our integration stepwidth 
is relatively coarse with = o.l. 

The most interesting regime of the ju 2 ] -curves is found near 
x •-> o and for small y (see Fig. 10). Indeed, as anticipated in 
section 2.2., the gradient of the computed curve becomes very 
steep there. If one extrapolated |u 2 | at x * o and y - o from 
data taken at y *f> o one would underestimate the actual value 
there. 

4. The relative importance of the edge region 

In the preceding section, it had been mentioned, that under 
almost all conceivable circumstances, a parabolic pressure field 
is created close to the plate edge. The relevar-je of this region 
to the shear layer excitation will be shown with a simple model 
(see Fig. 11 ) 



Fig. 11 Shear layer excitation by a monopole source. 
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It will be shown, what the shear layer velocity v 1 is far down- 
stream of the edge, if a pulsating (two-dimensional) monopolo 
source is located at different positions (a) close to the shear 
layer, (b) above the edge or (c) upstream of the edge (see 
Pig. 11). For the following calculation wo will only consider 
the amplified instability wave constituent and not the decaying 
one, because this is irrelevant at great distances x »y o /ui t Wo 
have from eg. (23) 


\ 


t , 


(90) 


u„ 


X » U Q /td 


A -i , 

e 1 


e" X 1 x v if dx 


with 


u 


(i+1) . 


If the upper boundary of the integral in eq. (23) is set equal 
to infinity, like in eq. (90), the total influence of the exci- 
tation is included. This expansion will provide the magnitude 
of the instability wave downstream of the interaction region 
with the monopole field, in section 2, we had already given the 
induced field of a monopole field near 5 semi-infinite plate. 

We have for v^ f (see also [12]): 


(26) 


If 


- /Hi— 7 

4ir q o 


Ax 


x+c o 

(x-x o ) 2 + y 0 2 


2 2 2 

with :r Q ® x Q + y Q and Q being the source strength Q » Q Q e 
The last term of eq. (26) can be split into two parts 


-iwt 


(91) 


x+r^ 


(x-x o ) 2 -y 0 2 


x-z 


x-z 


o - 1 


with 

(92) 


= x +iy *, = x -iy ) T - 1 +- 

n *♦ n ’ n n J r\ f 


i(r +x ) 


T = 1- 


1 ( r o +X o > 
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we are left with integrals of the type 
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(93) 


' . JL e U 1 • erfc /TTT , 

/x (X+T ) /t 
O 


The solution of the integral in eq. (93) can be found in tables 
for Laplace transforms or in the Handbook of Mathematical Func- 
tions (31# p. 302) by Abramovitz and Stegun. We end up with the 
following analytic solution *) 



but the coefficient governing the magnitude of these waves is 
fairly complex. We will# therefore, expand the solution for two 
typical cases, i.e., an excitation by a monopole source further 
away from the lip (-2— » 1) and an excitation directly at the 
lip of the semi-infinite plate. 

For the excitation at large distances, we have to expand the 
complex error functions for large arguments. For an arbitrary 
complex argument z we have for large z (31) 

(95) e z erfc z = — — M + 2 (-1 ) m - — - - - ~ — "] , 

/i?z *- m*1 (2z 2 ) m -J 

If we take only the first term of the series expansion eq. (95) 


*) This solution differs by the coefficient 1/2 from the solution 
given in (12), where this coefficient had been omitted erroneously. 
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wo find after some intermediate calculators P®®^ ^ 


(96) 


v 1 ■ 


Ota 

4U, 


‘o^S 


~ [' -ufc ( '- 2 72 - >] ■ 


This aquation can bo expressed also in terms of the distance r Q 
between source and plate edge and in terms of the angled (see 
Fig. 11). We find 


(97) 


v * — S_ sin * — ,/£ e x 1 x fl+ (4 cos 2$+ cos^)], 

1 2 *<^ "o n ' L ^ 1 


if we are interested in the modulus of v-j only, we have 


(98) 


ex q 

' v i'“ i^= sl ”t -i If’ f °*° ['* ^ “ s 2 »>] 

o * o 

wl* 

for » 1 and x » r , 

U ° 

o 


The first part of this equation resembles very much the excita- 
tion velocity in the neighbourhood of the plate edge 


(99) v, £ 


0 

2ir/r^ 



_L 

fk 


x « r 


One would obtain the following expression 

mx 

1100 ) |v,| . -2= sin • &- . °° 

2 "^ K R 

for {££ >> 1) 
vo 

if one calculated the excitation by the Darabolic field at 
the plate edge (eq, (99)) alone. Therefore, the expression in 
brackets reflects the additional interaction with the source 
field further downstream of the plate edge. The deviations from 
the pure interaction at the plate edge become small if er 0 /U 0 > 1 . 
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The conclusion from this is, that the parabolic field at the 
plate edge dominates if the source is further away from the 
edge. Now consider the situation shown in Fig. 11. In which lo- 
cation (a) , (b) or (c) of the source the interaction will it be 
the strongest? Equation (98) will give a clear answers at (c) , 
upstream of the shear layer! 

In that context an interesting question iB, how far an "exteriour 
excitation" can come from the turbulent shear layer downstream 
of the plate edge in a real flow situation, in our model, the 
pressure sources of the shear layer motion lie in the y - o 
plane. Therefore, in this model, no feedback from the downstream 
perturbations is possible (sin ^ = o) . However, in a real 
situation V 2 p = o is still valid outside the shear layer. The 
pressure sources are in a region of small O'. Therefore a very 
weak feedback of the downstream turbulent flow is possible. This 
consideration is not that naive as it seems at the first glance, 

0 ATI 

because the equation p = -2p is also valid in three 

dimensions. The source term on the right hand side might look 
slightly different in a nonlinear flow situation, but the con- 
cept of having linearly superposable pressure sources in the 
shear layer will not break down, because the pressure is a 
linear quantity in all our equations, and deviations of this 
linearity will occur only if the pressure perturbation is of the 
same order as the ambient gas pressure. Anyway, eq. (98) shows 
clearly, why shear layers are highly sensitive to perturbations 
(such as sound) coming from upstream and not very sensitive to 
perturbations having their origin downstream of the edge in the 
shear layer. 

Those who know recent experiments on excited jets know that an 
excitation close to the lip is very efficient. The preceding 
calculations did not consider this case, because it was assumed, 
that wr o /U o » 1. Clearly, this assumption does not include the 
excitation directly at the plate edge. On the Other hand, 
equation (94) can be also expanded for wr Q /U 0 << t, which would 
include the lip excitation case. To do this, we use the follow- 
ing expansion of erfc z for small arguments of z (see [31]): 
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( 101 ) 


erfc z * 1 - 


2 „-z 2 : z 2 "* 1 

“ !L 1 * 3 . ~ , (in.l) ’ 

n=o 


/it 


For small z, as in our case, we use only the first term of this 
expansion. It is no difficult problem, to extend the expansion 
further, but it will not provide much more insight. During the 
calculation it is also useful, to replace x Q by -x to avoid 
ambiguity problems with complex roots. At the end, X Q is used in 
the previous sense again. After some intermediate calculations 
we find 


( 102 ) 


v. - 2iL_ ^1 X 
2 U q 




oM cos(y oV - 2l p27 


sin f 


with X. » ~~ (i+1) and valid for ur o/U « 1 and “ X /U » 1 . 
1 o a a 

o 


For very small -+0 we find a very simple result 


(103) 


v, = e X 1 x 

1 2U„ 


This shows, that the source assumes its highest efficiency, if 
it is close to the lip, just above the shear layer. If we compare 
equations (100) and (103) in their relative efficiency, we find 
a coefficient n 


(104) 


sin I 

/#'• V 2 



This shows, however, that the excitation via the interaction 
from upstream (sin ^ = 1) and at a distance not too far away 
is not much worse than the excitation just at the lip. 

There is some interesting physics hidden in eq. (102). Assume 
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j y 0 \ 1 1 »+ o and i ! r Q ^ sin f >o. Then we have a source just above 
the shear layer at small, say, positive H Q , The source acts then 
ns a 5- function with the strength Q/2 on the shear layer. Q/2 is 
just the flux which penetrates through the shear layer plane, 
With this in mind we reconsider the general solution for the 
shear layer motion, eg. (23) 

X X 

(23) v, * - JL e X 1* [ e" x ix v^d* + A o X 2* J o” x 2* v 1f dx . 


For v.j£ m have hare 

(IQ5) v if * “ § 5 (x-x Q ) at y * +o 

with oq, (23) wo have the complete solution at once 

[e x l x • q“ x 1 x Q - a X 2 x . e ~ v 2 x o] 


UQG) 


f, ■ R ,3pW 

QuV 


V.J * 


2U„ 




If both influences, the "local" #- function influence and another 
term of the parabolic typo are in competition, one sees, that 
the "local" influence can become weak in comparison to the para- 
bolic typo influence. In particular, if the source is further 
away from the plate edge, the parabolic influence at the edge 
has piled up to such high magnitudes by the exponential amplifi- 
cation so that a local influence (see position (a) in Pig. VI) 
would only contribute a term of the order 

I' 07 ' o'l«-V . 

For the pure "local" excitation with I A 1 x (J | -*■ o we have, on the 
other hand, a simple analytic solution which might be utilized 
if the free shear layer is really excited ‘just at the plate edge 


') with G 1 » G 2 « o 
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Also the u-disbributlon in the whole ambient field can be written 
in closed form, it contains exponential integrals of complex 
argument for the induced field of the instability waves truncated 
at x * q (v 1 « o for x <o) and of the induced field of the source 
The details can be seen in appendix C*) whore just this calcula- 
tion is given in another context. For distances further down- 
stream, we have a ultra- simple equation for the magnitude of the 
instability waves, where jv^| » 1 v 2 1 » |u^ | * [ u 2 | . We find 


(109) 



at y * -o 


with the usual decay in t y direction 


( 110 ) 



x » u Q /u 
y '+ o . 


2 

Q is the volume flux (say, in m /a) of the excitation source. In 
a real situation with an arrangement like the one shown below 
(see Fig. 12 ) we suspect, that more than half of the volume flux 
Q penetrates through the y * o plane. Therefore, the efficiency 
might be even slightly higher than suggested by eq. (110). 


*) For numerical computations in the potential field outside the 
shear layer, it is also recommended, to use the calculations in 
(36). 
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Fig. 12 real edge excitation configuration. 


5. Two streams with different densities 


As we will see, there is no simple transformation to convert the 
results of the case with a single stream (U^ * o for y > oj 
U 2 = U G for y < o» P 1 - P 2 ) to the case with two streams 
02+0} p 1 *|»p 2 )‘ Nevertheless, it will turn out tfrf/t the results 
of this section will be still very similar to those of section 2 
However, the mathematics is much more tedious. Since it would 
not provide any specific insight to show all details of the 
calculations, we will provide here only the basic ideas and 
steps. 


We start, as in section 2, with the condition Of equal displace- 
ment h at both sides of the shear layer 
f 


(111) < 


V s ill + n 
v i at + °i ax 


v _ 111 + u 3h 
v 2 at u 2 ax 


J . 


_ j , . 4, 

By using h'e’ and by eliminating h from equation (111) we 
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«UR 


obtain one condition relating v 1 and v 2 


( 112 ) 


U„ Ov. U i »V„ 

v i * * i -sr ■ v 2 * 1 i tnr 
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This equation holds for all x and is also valid for "j 4 p 2 # 

since it ia a merely kinematic condition. To create a second 

condition for and v^ we roly again on a consideration of the 

pressure field, since the density jump is assumed to occur in 

the shear layer as well as the mean velocity jump# we obtain 
2 

again v p » o outside the shear layer and with exception of the 
location of tho oxter, tour source. Therefore# we can again use 
the second Euler equation (2) 


_ av.. apt 
13 1 ( “ iwv i * u i ar> P “ oT 


(113) \ 


n 


°2 < " Uv 2 + u 2 1T ] 


ay 


The procedure to split tho pressure field into two constituents# 
the forcing terms (index M f") and the shear layer terms (index 
"a") is still possible. The "£" pressure gradient is continuous 
through the shear layer and the n s" pressure gradient ia sym- 
metric with respect to tho shear layer. We have 


(1 14) 


p i ( “ iwv if * " fl 2 ( — d. is v 2 £ + u 2iT* > 


Ov. 3V 5 „ 

»1 ♦ “l-TT 1 ■ -l>2 | - 1 “'' 2a 4 U 2-#> 


whore v^ « v^ v^ 3 ; v 2 * v 2 £ + v 2g . Equation 03) can be con- 
densed into 


. U, Dv, U, Dv 0 v 

(11D) (v-, * -1* ax 1 ) + v 2 + ir T o^T " 2(v 2f + i V "ilr^ * 


U 2 3v 2f 
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Since it: is more convenient to calculate the forcing pressure 
field (which is still parabolic around the semi-infinite plate) , 
we prefer to express the right hand side of eq. (115) in terms 
of the pressure fipid. We have from the second Euler equation 


(116) 


U„ DV 


2f 


+ i 


2f 


9X 


1 3p 2f 
p 2 ui Dy 


where the pressure gradient in y-direction is taken at y = o, 
We will abbreviate in the following calculations 


(117) 


P * 


21 0p 2f 
P 2 « ^ 


o 


) 6 



Consequently, we can qrite eq. (94) 


(118) 



with the right hand side being known and proportional to 1//x. 
Equations (118) and (112) can be used to provide a single non- 
homogeneous differential equation for v^ 


(119) 


21 Ur 

(1 + B) V 1 +---^ 

! wJ 


(1+oB) 


^1 

dx 


U 2 „ d 2 v. ioU 
_! (l +o 2 0) -/=£+— 
in dx 


2 d£ 
dx 


The homogeneous solutions (Helmholtz solutions) are 
(120) f 1 = C 2 eH x ; (f> 2 * C 2 e X 2* 


with 


( 121 ) 


U 2 1+o 3 


i ( 1 +0 0) i ( 1-o) /0 
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(122) 4§ * |o x 1 x | o" v l x Uiidx~o. x 2 x J 0 “ v ?* x RSdx 


whom "rs" stands foe the right hand side of eg . (113). With the 
same arguments as those given in appendix 2 for r - o and n ° 1 
wo can sot e G 2 a o» We can further evniuate oq. (122) 


(123) v 


1 2U 2 /d thA* 


(ii-io if)e'' ,x | o" v ' x pdx “(i-io/R)e x 2 x o“' 2> Vxlx 


for a parabolic pressure field around the semi-infinite plate wo 
have from oq. (33) in section 2.1, 


(124) 


21 , ;lp 2C 
v'y 


21 , Ap 1 2 


wp 


2 4 tlx 


Consequently, wo can again relate all quantities in our flow 
field to the excitation, characterized by the pressure difference 
Ap .| 2 at both sides of the plate at. the distance 1 from l-he plate 
edge. Therefore, we can introduce dimensionless quantities simi- 
lar to those used in section 2: 


(125) [ v 


V, a 


o V p 2 > ' u,D 2 1 
A[) i 2 ' a A P;i2 


v j • pgt'VDgi 


v, 2 « “V 

1 < A 1 t-o i 

- |i(1 + oi<) 

y„ u)X . 

X « , O 

u 

a , U es 

G 2 

U 2 


n , 


in the expression Cor \ the positive sign in the brackets cor- 
responds to \ 1 , the negative sign to respectively. With 
those abbreviations wo obtain the final result for , after 
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#1 w. » 1 v , « „ v A V H j 1 *1; 1, y |iv . \ 1 X..,, p ,|i r* 

1 1 « u 3 * v * ut ' rw - - * £x * j 1 i7< • ?■ 1 1 1 «j \ * x ~ 

\ 1 4 *iUU**iO 1 >•? 

] (X S Of V * H)J 






To obtain u 2 WP of ' n either derive a similar equation an eg, (119) 
from the Home initial, equation* dial ant) (118) o* utilise eqs, 
( 112 ) ami (118) toy ether with eg. (124) to obtain \*> from Vj. 

We end up with 


{ 1 * U ' 1 Vi ”o r f V \ ,, se 1 zMl&o K 2 *orf/\ 2 
1 

■ (x > o, y » -n) 


(127) * v., - - •» ' ^ - , s v , . 

4(1 't v»“8) ' > \ 


He* itto.tt the differences between v 1 and v ? which resemble those 
in section 2, we have a different coefficient (»8) in the deno- 
minator in front of t he bracket a. Thin leads to different vol o- 

ciU.es of the instability waves on both sides of the shear layer. 
Nevertheless , the ill sphu'ement tv at both aides of the shear 
layer in live same, as postulated in the beginning, 

Before we determine the corresponding u-diatribations we will 
discuss the behaviour of v ( and v 9 for amall and large x, For 

small k the expansion of eg, ( 1 2 6 ) and (127) yield 

U -* ** . JVX 

2 Hv*s 


'^1 5 

in both canon the envelope of the displacement h la h « x , 
consequently the mean flow leaves the trailing edge of the plate 
tangentially, which is equivalent to the Kvitta condition . For 
d.| *o and - a *o higher order terms in the expansion' for v ( , 
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,<-/x 


(x * o) 






, * 




•Cl 


"i * * 


take over, ao that h is ati.U h « x 1 - 5 , ORIGINAL PAQI IS 

OF POOR QUALITY 

For large x wo obtain for v., and v 2 

i»-ne^ > x 


(129) v 1 «» 


4v H (1-inv‘B) i \ , 


4»> 


(X + «) 
(130) v 2 = - 

(X ► '«) . 


VTtO 


\ ) X 


4 (1“ivtvl<) > \i 4>x 


2 

1 +f 


X 
1 + 1< 


l * 


Wo again have a Hpliriholta-typo solution (first tonus of egs, (129) 
and (130) ) in tho far field and another contribution due to the 
excitation being proportional to 1 / < x . The interesting thing is, 
that nature takes here an av erage value of the densities repre- 
sented by the coefficient 2/1 +p. However , the Helmholtz- typo solu- 
tions are different in their magnitude by the coefficient 
in the different flow regimes. 

The u-velQCitios are calculated in exactly the same way as the 
V* velocities above, with an approach identical to that used in 
section 1, we. end up with 


(131) 


(132) 


(133) 


(134) 


u ^ * + 


4>1i(1+n B) 


1 ->■ .1. o if \ ,x _ 1-ln/i< \£X 


* 


(x > q) y ® +o) 

U*_ ■ 

2 4 ( 1 + o 2 B ) t_ v'\~ 

(x > o,* y * -o) 

/** 

* V /’II 


^2 


Q 0 M X + c : '2*”| 

*“ A-) * \ j 


1 


4. / B(1+y B) 


lli X B e ^l x erfc ,^x- JXXi c ''2 x et . fc ^1 

- sq >q; 


(x < o; y * +o) 


a - i-4- (l±i£iiS d X 1 x arfc 2 * er f c STS 

2 4(1+0 e) sq iT£ 

(x < o> y = -o) 
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For large x only the amplified Helmholta-solutions survive (terms 
with X 1 in eqs ( 1 3 1 > and (133) ) .The modulus of th.-ui amplified 
waves is for x » 1 


(135) 

i r. , _ 

-o) X 

1 4/| < 7(1 + 0) (1+0*3) 

x » 1 

and 



(136) 

iu 2 |= /s • ! u n 1 



x » 1 



If we return to dimensional quantities, we have for eq. (135) 



6, Theoretical considerations related to experiments 

Parallel to this theoretical investigation also experiments were 
carried out. On these experiments we will report later. During 
the experiments a number of typical problem-? arose which can be 
settled theoretically rather than experimentally. 

A point of particular interest is, that the phase speed of the 
instability waves in the experiments was lower than that of our 
theoretical model. Phase speed and wave number are inversely 
proportional to each other. Since the wave number is coupled to 
the exponential decay rate (of the fluctuations) perpendicular 
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to the shear layer, this has a significant influence on the in- 
duced field (see eg, (68) in section 2,2,), Conditions which can 
change the phase speed are : "overshoot" of the mean velocity 
profile (section 6.1,), finite shear layer thickness and entrain 
ment effects (section 6.2.), Finally, we have to consider the 
effects of additional walls in a test facility. We will consider 
the case of a shear layer in the symmetry plane of a rectangular 
channel (section 6.3.). 

6.1. "Overshoot" of the mean velocity profile 



Fig. 13 Interaction between nozzle and screen. 

In Fie- 13 it can be seen how an "overshoot" of the mean veloci- 
ty profile can be (inadvertedly ) generated. A flow penetrates 
through a screen of high resistance at (almost) constant speed. 
The following contraction in a nozzle leads to different accele- 
ration of the different streamlines. In particular, the regions 
closer to the wall are accelerated more than the flow in the 
center of the nozzle. This creates the "overshoot" of the mean 
velocity profile. Since it occurs often in experiments it is 
worth investigating its influence with a simple theoretical 




model. As we will see, the deviations due to this effect cannot 
bo neglected, , 
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The model which we are going to consider can. be seen in Fig. 14 . 
Since we cannot use our pressure gradient considerations as in 
the preceding section we will proceed with a slightly different 
approach. In addition, we will restrict ourselves to the dis- 
cussion of conventional spatial instability waves extended from 
x s to x a +m, The aim is to find how the wave number and the 
amplification rate will change if we introduce a mean flow gra- 
dient (see Fig, 14) . 

We start with the same basic equations as in section 2., i.e. , 
the linearized Euler equations (1) and (2) and the continuity 
equation (3). If we subtract the x-derivative of (2) from the 
y-derivative of (1) we can eliminate the pressure and obtain 


(138) 


. .au av . 
' i<J ( ay " to* 


+ U 


ax 


au 

ay 


— ) 
nx' 


+ V' 


a 2 u 

ay 2 


This is what remains from the inviscid vorticity transport 
equation if we consider a linearized fluid motion which is 

”' i ' lllt ) and takes place in a parallel mean 


harmonic in time (« e 
flow with velocity U 


(y) ’ 


If we assume a potential flow with 


u ~ a^/ax and v = &$/ay we have from eq. (138) 
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This means that we can use a potential flow approach in regions 
wherea 2 U/Dy 2 » o. This condition is fulfilled outside the shear 
layer. Ah both sides of the shear layer wo have to match the 
displacement h and the pressure p* as before. We use the follow- 
ing ansatz for the stability waves 


for y > o: 0^ * o ; ^ 


Ae ia (x+iy) 


(140) < 


dU 


for y < U 2 * U Q * gy • y ? * a 


Q^ia (x— iy) 


Since $ is a function of x + iy»the continuity equation written in 

2 

terms of the potential, i.e., v <j> * o is satisfied. The potentials 
fulfill also the boundary conditions at y *'<*>, because the in- 
duced velocities become zero there, a is the complex wave number 
defined as usual in stability theory. A and B arc magnitudes of 
tho potentials which are still arbitrary at this stage of the 
calculations. 


The kinematic condition of equal displacement h at both sides of 
the shear layer is, as in section 2: 


(141) 


V 2 


* V^ + i 


u> 3X 


Using eg. (140) we obtain 
aU 

(142) B « -A (1-~) . 

w 

The second condition is, that the pressures or the pressure gra' 
dients should be equal at both sides of the shear layer. For 
y > o we haye from the first Euler equation (V): 

ap. 

(! 43) = i P uiu 1 
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and for y < o we obtain from the same equation (1) 


(144) 


3p, 

rt- 


-p (-iuiu 2 + U 


3Uj 

o alT 


dU 

3y 


V 2> 


If equations (143) and (144) are set equal and with eq. (140) we 
have 


( 145 ) 


A * B (1- 


V 



Equations (142) and (145) provide a quadratic equation for a with 
the solutions 



if we compare this to our previous nomenclature, using 2 =ifl i 2 
we find 

i’ 47 ' »i,2 [<■ M t\/T- <#*7" . 

In comparison to the situation, which we had before with 
dU/dy - o, we have now an increased wave number (imaginary part 
of eq. (147) ) , which will cause a more rapid decay of the induced 
field in the y-direction. The deviations can be, for a typical 
situation in an experiment, in the order of 10-20 %. On the other 
hand, the influence On the amplification rate (real part of eq. 
(147)) is much weaker and will be in the order of typically 
1 - 2 %. 


6.2. Entrainment and finite shear layer thickness effects 


If one considers the effect of entrainment alone for an infini- 
tesimally thin shear layer, one can extract the relevant informa- 
tions from section 3, equation (121). For constant density, the 
wave number is 


(148) 


Im (X) = ~ 


1 * V°2 

0 2 1 -mu , /u 2 , 2 


b5 


, jluu 
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For small g « U^/U 2# say cr « o.l, we have a wave number increased 
by 5 9 % . This enhances also the exponential decay rate in the 

y-direction by 9 %. 

As expected, the downstream growth rate decreases with increasing 
g. We have for the growth rate 

<149) Re <\) - * ■ 1— ? — y . 

U 2 1 + (U^Ug-r 

For, say, a = o.l we have a growth rate decreased by # 11 %, 

These considerations are valid for a "thin" shear layer. 

We cannot really neglect, however, the effect of finite shear 
layer thickness. A typical parameter of the stability analysis 
of shear layers of finite thickness is the Strouhal number 
S e B 0 *£/U 2 where f is the frequency and Q is the momentum thicks 
ness. The momentum thickness is defined as 

(150) 0 = f JL (i- JL)dy , 

°2 °2 

For a typical shear layer profile, which might be approximated 
with good accuracy by a tanh-profile we have the situation dis- 
played in Fig. 15 



Fig. 15 Shear layer profile 
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The amplification rate and the phase speed, according to Michal- 
ke [6 } are given in Fig. 16 

As we see from Fig. 16, the deviations in the wave number are 
much more significant for increasing S 0 (and thus deviating from 
the "thin" shear approach) than for the amplification rate. 
Therefore, the finite thickness causes an increasing wave number 
and enhances again the decay of the induced field in the y-di- 
recfcion. This influence can be quantified with the numerical 
data available in the literature (6, 37). Also the combined in- 
fluence of entrainment and finite, shear layer can be extracted 
from (37), where calculations are given for a shear layer of 
finite thickness between two streams of different velocity U 1 
and Dj. 

6.3. Channel effects 

In an experiment, one cannot offer a facility with infinitely 
extended streams at both sides of a shear layer. Therefore, we 
will have to consider effects due to the finite dimensions of 
the two streams. The experimental setup which can be modeled 
with the least complication is that of a free shear layer in the 
symmetry line of a two-dimensional channel. At the rigid walls 
the velocities and pressure gradients normal to the surface have 
to be zero. Thus, the boundary conditions at both sides of the 
shear layer are the same. Consequently, our basic approach which 
takes advantage of splitting the pressure field into a symmetric 
and an antisymmetric part, is still valid. This has been shown 
already in a different, more involved way in [12], Therefore, 
the general solutions for the v-velocity component are still 
valid (eqs • (23) and (123) J . However, we will have to discuss the 
deviations in the actual velocities and pressures which will 
occur in 

(i) the excitation field (v-velocities) 

(ii) the resulting shear layer motion (v-velocities) 

(iii) the induced u-velocities 
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There is no doubt, that the structure of the excitation field 
will be still parabolic if we are close enough to the edge of 
the semi-infinite plate, Therefore, there will be particular 
conditions when the solutions calculated in the preceding section 
are still valid. The restrictions of these conditions will be 
discussed below. We will calculate the fluctuating flow (without 
shear layer and mean flow) in a simplified test section, which 
can be seen in Fig. 17 , left hand side. 


2 a • plant Z b * pioni 



z a -plare (physical) z^-plane (image) 

Fig. 17 Simplified test section and its conformal image. 


The test section is simplified in so far as the nozzle c. I the 
settling chamber on the left hand side have not been considered. 
The excitation is assumed to be generated at great distance on 
the left hand side inside the bisected channel by two sources of 
opposite sign**. This would be a good model for the situation 


** A similar procedure as the one described here could be used 
if the source were in the lower and upper wall, created physical- 
ly by vibrating plates. 
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with sound excitation by loudspeakers connected to two settling 
chambers upstream of the bisected channel. 

If the sound wavelength is much greater (say, at least ten times) 
than a typical dimension of the problem, (e.g., the channel width), 
one can describe the alternating flow in the vicinity of the 
plate end by incompressible flow theory. In that case, we can 
use the very powerful tool of conformal mapping to solve the 
problem completely. For configurations like the one shown in 
Fig. 17, one can use the Schwartz-Christof fel-transformation to 
map the flow in the interiour of the channel into the flow in a 
half-plane (Fig. 17, right hand side). Our channel and splitter 
plate walls will be stretched into the straight line y b = o 
(see Fig. 17). The open ends at the left hand side of the bi- 
sected channel will be mapped into the two points at x b ■ +1 and 
x b = -1 on the x^-axis. If we install a positive source at x b =1 
and a negative source at x b = -1 we have already the image of the 
flow field in the z^-plane. 

The mapping function relating the z & and z fa planes can be found 
in Kober's dictionary of conformal representations [38], We have 

z b - 

with z a = x a + iy a (physical plane) and z fa = x b + iy b (image 
plane). This mapping function corresponds to the specific dimen- 
sions given in Fig, 17. We introduce the complex potential 
(|j = (j> + i¥ where the potential function <(> and the stream function 
are defined as usual in plane flow problems. The complex potential 
of a source at the origin is = In z. Thus we have for the 
potential of the positive and negative source in the z^-plane, 
located at x b = ±1 

(152) [J] = A [in (z b -1) - In (z b +1)J 


where A is a constant (proportional to the strength of the 
sources). With the transformation equation (151) we have in the 
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053) (Jj = A [In ( Vlte^a-i) _ i n < -/we^a+i) 


We are first interested in the v-velocity (in our previous nomen- 
clature v^£> at the center line of the channel. The complex velo- 
city is 

054) I®- « u - iv . 

az a 

Taking the derivative of eq. 053) we obtain 

(155) u - iv = — A — -■ . 

l/l+e Za 

The center line of the channel is located at y„ * ir; there we 

. a 

have z = x„ + iir. With e lir = -1 we obtain 
a a 

(156) u - iv = — -- — at the center line (C.L.) , 

Vl-e Xa 

For x a < o this is purely real, which means that y = o there. 

For x a + we have u c L + u o , the homogeneous velocity at 
large negative distances in the bisec-ted channel. Therefore is 
A = u^. Equation (156) becomes purely imaginary (and u = o at 
the center line) for x_ > o. 

a 

Before we write down our equations in real physical quantities 
such as u m , we have to relate x also to the only real physical 
dimension of the problem, i.e., the channel half-width d (see 
Fig. 17). Therefore, we will replace x a by a dimensionless quan- 
tity, i.e., The coefficient ir occurs, because in the initial 

coordinate system the channel width had this value. We end up 
with -.he following equations for the velocities at the center 
line y = ±o: 
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For small positive x we have 

1159) y - u. • VI ' 'tji ■ 

1 

This is the expected /fa behaviour of v. On the other hand, it 
can be seen from eqs, (157), that the deviations from this 
behaviour can become significant. 
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Fig. 1 8 shows a log-log plot of eq. (157) as well as some own 
experimental data*). As we see, the theory works quite well. A 
second question of major importance is how far the parabolic 
pressure field is valid upstream of the plate edge (see Fig. 19 ) 
In particular, wo are interested in Ap^/Zl as a reference 
quantity. 


rmrrn 
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/• / j y—jP -const 
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Fig. 19 Reference pressure locations in the test section 

Since we can expect a parabolic pressure field only in the near 
vicinity of the plate edge, we would run into difficulties with 
measurements of Ap^/ZI for the following reasons: The pressure 
difference of the excitation field Ap^ becomes small at small 
distances 1 from the edge. This leads to poor measurement 
accuracy. Close to the lip the signals from the (symmetric) 
pressure field induced by the shear layer dominate and determine 
the setting of the measurement system. However, for the measure- 
ment of the relatively low magnitudes of the pressure difference 
the measurement system is then poorly adjusted and the measure- 
ments become unsatisfactory. The situation would improve con- 


*) The points correspond to pressure gradient measurements 
(3p/3y) taken with a probe microphone of =1.5 mm resolution, in 
air at room temperature (22°C) and at a frequency of 200 Hzj 
d was 75 mm. 




siderably, if we could take data of ip (still at both sides of 
the splitter plate) further upstream where the pressure levels 
of the excitation are higher and the induced (symmetric) pressure 
level of the shear layer would be lower. Consequently, it would 
be desirable to determine ap i 2 //I# which refers to the parabolic 
field near the lip, out of a measurement &p further upstream. We 
will achieve this with the same method of conformal mapping which 
we have used before. As we have discussed in section 2.1., the 
pressure is proportional to the potential for a fluctuating flow 
field with no mean flow. In analogy to eq. (153) we can set 


(160) p = B-Re •{ In ( V”+e z a- 1 ) - ln( VT+e^a+l)} 

Where B is a constant which we will determine later. We restrict 
ourselves to the center line of the channel with y a = 11 and we 
replace later x by Tix/d, as before. We find 


VT Trx/d . 

(161) p = B • Re( In (-ilZ§ Z_L) Y 

Vl-e^+l ' 


The denominator of the argument of the logarithm will become 
negative. Since In (-1) = +in. this negative sign will only con- 
tribute to the imaginary part of p which we do not consider here. 
We have therefore as the real part of p 

t ,/i "3x/'d 

(162) p = B • In (i-^1-2 ) . 

1+ Vi -e' rX/d 

For small x we obtain 

(163) p = -2B /77d • /^x 
x ■+ o . 

This is the expected parabolic pressure field near the plate 
edge (see also eq. (33) in section 2.1.). On the other hand we 
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-LI B lE_ , 
/I /=H 


because the pressure difference Ap at both sides of the semi-in- 
finite plate has double the value of the pressure at one side 
alone. Using equations (163) and (164), we can determine the 
constant B: 


( 165 ) 


B = 




Ap 


12 

rr 


Using Ap = 2p , equations (162), (165) and replacing -x by U (see 
Pig. 19) we obtain as a final result 


( 166 ) 



This provides the desired relation between the pressure difference 
Ap at both sides of the splitter plate, measured at an arbitrary 
distance L from the plate end, the channel width d and the 
reference quantity Ap^/^^ (see also Fig. 19), Eq. (166) can be 
expanded also for L >> d which gives 


(167) 


AP 


12 


A£ 

L 


/I 

at (L » d) 


/d 


_2_ 

/if 


One should keep in mind, however, that L should be much smaller 
than the acoustic wavelength, which restricts the applicability 
of eq. (167). It is, therefore, more likely that eq. (166) will 
be applied in an experiment. 
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6.3.2, Effects on the resulting shear layer motion 

We restrict our considerations to the single stream case where 
we have a mean flow U Q in the lower half of our channel. The 
general solution for v is still valid, as mentioned above 

e “*2 x Vl f dx 


with X. 0 * — ( i * 1 ) . For the vertical excitation velocity v 1f 

' t * n 1 x 

we have to use now eq. (157) 


(23) 


V)- 


gj 


, x i x 


c^^v^dx + 


o X 2 x 



(157) v 1f 


We will focus our interest on simple results downstream of the 
interaction region. In this situation we can neglect the second 
term of our general solution (eq. (23)), because it represents 
a rapidly decaying wave. As we know from our previous results, 
the interaction region is limited to a distance of less than one 
wavelength of the instability waves. Since v 1f is now decaying 
even more rapidly, we can replace the upper limit of the integral 
in eq. (23) by infinity. This will provide completely reliable 
results for — » 1. So we are left with 



(168) 


v. = -u 

1 o 


G X 1 X 


■X-|X 


dx 


for (— » 1), 
°o 


/TtX 

°/e d -1 
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The integral in this equation can be rewritten in the following 
form 


h* , f G ’ 1 2d' 

. dx = ~ dx • 

^x J -Tb< i 

d_ 1 o ( 1-e d ) a 


(169) 


/e 


This latter integral can be found between the Laplace transforms 
in [39] . We obtain 


(170) 


1. = 


1 X 1 d 

Bill * -f) 


>. 


where B (x,y) is the Beta function (or the Eulerian Integral of 

the first kind), defined as 
i 

(171) B(x,y) - | t x_1 (1-t) y_1 dt ■ 

o 


The function B(x,y) is related to the f function 

(172) B (x,y) = B(y,x) = r (x) ‘ r (y) 

r (x+y) 

The P function (or the Eulerian Integral of the second kind) is 
defined as 


(173) 


<z) 


*/ 


e _t • t z-1 dt 


with Re (z) > o , 


Using eqs. (170) , (172) and r (-j) » /if we find 


1 1 1 d 

r (A + — !— ) 

(174) vi = -u • — • — • --2 - ... I 

" U n /ir X.d 

° r (i + — !— ) 


Before we proceed with the discussion the r functions, we have 
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to establish a relation between and our reference quantity 
A p^ 2 /^l* We have to expand eq. (157) for small x to obtain the 
parabolic flow region of the excitation, to which the reference 
quantity Ap^/*^ belongs. We find 


< 175 > v if - M* 


on the other hand we had (eq. (36)) 


(176) v 


if ~ " fk 4pu) ’ /i 


Comparing (175) and (176), we can replace u in (174) 


1 X l d 

,,771 v, - • A • IS ■ • - T-* 

' /o ° 

It 


or, if we use a dimensionless v 1 , as in section 2 


V 1 = 4 


r a + -U) 

l 

X.d 

r < i + — H 


X,x 

* e 1 


This is a new analytic solution for the excited instability 
waves in a channel at ^ » 1. The r function is tabulated for 
complex arguments in [3?], However, it is difficult to see what 
eq. (178) really means. In particular, we are interested to see 
what the errors are if we use the equations for a free shear 
layer without walls (section 2) in the present situation in a 
channel with walls. We want to see by a quantitative calculation 
under which conditions the previous approach will fail. It makes 
sense, therefore, to expand eq. (178) for large values of X^/u. 
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We will expect then, for decreasing values of A^d/ir (which is 
proportional to the channel width divided by the wavelength of 
the instability waves) an increasing deviation from the previous 
approach . 


Fractions of the r function can be expanded for large arguments 
in. the following way (31) 


(179) 


ftfrli ’ * a " b 0(z - 2 ,i 


Using this expansion, we end up with: 


(180) 


°1 = 1 


~.A . M- • ^2: p J '1 x 

4/T+T 8 0- + 1 • 


for — » 1 and — » 1 . 


With exception of the term in brackets, which quantifies the 
deviation caused by the channel walls, we have obtained the same 
result as in section 2 (eq, (49) , first term of the instability 
waves)/ We are also interested in the modulus of v^ , which is 


(181 ) 



/tt 

4-V2 


< 1 - 


JL V 

16 ’ ed 


for u)x/U Q » 1 and ud/U o » 1 , as before. The abbreviations in 
eq. (181) are the same as in section 2, i.e., x = wx/D and 
v^ = v^ pY«u o l/Ap^ 2 • Equation (181) is simple enough for a 
straightforward error analysis. Consider a simple example: 
d = 75 mm, w = 2u • 82 Hz, U Q = 12 m/s. We end up with a value 
for | v 1 j being about 6 % lower than predicted for a free shear 
layer without channel walls. This deviation reflects the in- 
fluence of the limited size of the parabolic excitation region 
around the end of the semi-infinite plate. 
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6.3.3. The induced u~velocibles OF POOR QUALITY 

It would be quite complicated to calculate the whole u-velocity 
field, Therefore, we will restrict ourselves to “ » 1 where we 
have to consider basically deviations from the U o field computed 
in section 2. It is quite obvious, that close to the edge of the 
semi- infinite plate the flow field Will be the same as in 
section 2, If we consider the field upstream of the edge, it is 
reasonable to assume that the field is dominated by the excita- 
tion. Consequently, at sufficient distance upstream of the edge 
we will obtain u * ±u m at both sides of the plate, In the inter- 
mediate range upstream of the edge, but closer to the edge, the 
field becomes more complicated. For the upper region with no 
flow, the u-velocity field can be derived from the known pressure 
field in a straight forward way, using the first Euler equation. 
With other words, u is proportional to the x-derivative of the 
pressure. In the lower half of the flow field, the first Euler 
equation produces a first order differential equation for u, 
which has to be solved. The same situation arises, by the way, 
also for the parabolic flow field in section 2, with the differ- 
ence, that the resulting functions (error functions) are there 
easier to deal with, 

The situation far downstream from the plate edge becomes much 
simpler. Here, the amplified instability wave dominates the ve- 
locity field. The effect of the channel walls can be taken into 
account by a simple imaging technique (see Fig, 20 ) 



iig. 20 Imaging technique to model the channel wall. 
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The decay of the induced velocities in the -y-direction is pro- 
portional to a^hdii. We Know then, that the u,-velocity at the 

II " m *)m4 

shear layer will o be increased by an amount (1 + e - ;r -) . Using 
the fact, that at ~ » 1 we have | v 1 | * |v 2 | * | u 1 1 « u o |u 2 | for 
a shear layer U ° without walls, wo obtain for the situation 
in the channel 


2wd 


(182) 


I V 1 


(1+e u ° ) 


(^ » 1 ) . 

5 o 


To give a numerical example, we take the same data as before, 
namely i d B 75 mm, w * 2n • 82 Hz, U Q * 12 m/s. We find a very 
weak influence of the wall of 0.16 %. This influence will become 
even weaker because the decay rate in y-direction will become 
higher in a real experimental situation, as discussed in 
sections 6.1. and 6.2. We end up with the conclusion in this 
section, that the enhancing influence of direct reflection on 
the u-velocity at the wall is usually much smaller than the more 
indirect decreasing influence via the change of the excitation 
field. 


7 . Speculations 

One may contemplate about the question on how far the present 
approach can be stretched to really include the excitation of 
a shear layer with small, but finite thickness. In section 5 
(two stream case) we have considered an example where the ampli- 
fication rate and the wave number can depend, albeit differently, 
on exteriour parameters like a = u.|/0 2 and (3 = p ^/P 2 * It is 
tempting to write an equation like eq. (119) where the coeffi- 
cients on the left hand side are made up from amplification 
rates and wave numbers obtained from numerical results of the 
stability theory of infinitely extended shear flows, such as 
[6] and [37] . It is not dear, whether or not this idea leads 
to a viable and correct analytic approach, However, it would be 
interesting to see how far the concept can be generalized, that 
the resulting excited instability wave is a finite Laplace trans- 
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form (such as eqs. (37)and ( 1 23) ) of an exteriour excitation 
field. In addition, extensions of such an approach to other 
configurations like jets or even to the excitation of instability 
waves on a flat plate might be conceivable. 

Another relevant question is, how far free stream turbulence 
(Or other types of turbulence) and its wavy predecessors which 
we investigate here, are "natural" or dependent on exteriour 
forcing. We might consider whether or not an influence from the 
downstream nonlinear development might provide a sufficient or 
even dominant feedback. Our pressure source approach is still 
useful to provide a feeling for this situation. The pressure 
linearization is the last to break down with increasing fluctua- 
tion magnitude and it works nicely for subsonic flows. The 
effect of the nonlinear flow evolution on the pressure sources 
is this: During the vortex roll-up we have a limitation to the 
magnitude Of the pressure sources. The reaction to the trailing 
edge is necessarily very small*) , because the pressure sources 

are still only located in the shear layer (outside the shear 

2 

layer we still have V p * o) , i.e. , in the symmetry plane and 
we get a situation close to creating a stagnation point of the 
induced flow at the trailing edge (see also section 4). The 
feedback is therefore strongly dependent upon the local displace- 
ments of the shear layer. Therefore, a feedback is likely to be 
a nonlinear effect, and it is very weak. Thus, it can be easily 
dominated by sound waves or convected vorticity, as shown in 
numerous experiments, Nevertheless, we cannot completely reject 
this mechanism, because the upstream influence decreases in a non- 
exponential manner in the upstream direction whereas the down- 
stream amplification of the instability waves is exponential . 
Under these circumstances (and with an unstable shear layer), 
there is no Way to keep the flow laminar even if there is only 
Brownean motion as an initial excitation somewhere. 

There is, however, another way in which feedback may come into 
Play. 

*) The situation changes dramatically If a rigid body, like an 
additional wedge, Is introduced into the flow field. Then, an 
asymmetric pressure field is created which has a strong influence 
on the trailing edge. 


All our previous considerations were based on symmetric boundary 
conditions. To provide a feeling for slightly asymmetric eon- 
ditions we consider the geometrically symmetric configuration of 
section 6, where we excited a shear layer in the center line of 
a channel. If the conditions upstream of the splitter plate are 
not symmetrical (caused, say, by a difference in the acoustic 
impedance in both legs of the bisected channel), a pressure 
sources in the shear layer creates an asymmetric pressure field 
in the channel and hence a parabolic local field at the edge of 
the splitter plate, In this way, in particular if there is 
resonance in one leg of the bisected channel, a linear feedback 
is conceivable. 

8. Conclusions and survey of essential results 

The acoustic excitation of a free shear layer shed from a semi- 
infinite plate has been studied. The shear layer is assumed to 
be infinitesimally thin. The present consideration is limited to 
a low Mach number mean flow and to the situation where all typi- 
cal dimensions of the problem are small compared to the acoustic 
wavelength. This leads to an incompressible flow approach 
relevant to the interaction region near the end of the plate. 

The basic contributions and findings of this paper are listed 
below . 

1) The basic shear layer - sound interaction model is out- 
lined. The solutions are obtained by splitting the pressure 
field into two constituents: 

(a) the forcing pressure field of the exteriour excitation 
and 

(b) the pressure field which is radiated by the shear layer 
itself. 

The consideration of the pressure gradients perpendicular 
to the shear layer leads to an equation for the v-velocity 
components v^ and v 2 at both sides of the shear layer. A 
second equation for both components is derived from the 


condition, that the displacement h on both sides of the 
shear layer should be equal. Both equations can be con- 
densed into one nonhomogeneous differential equation for 
v.| . A general solution for this equation is derived. For 
a parabolic excitation pressure field around the plate, 
which is the most common case, a special solution is given. 
So far the basic approach of this paper follows a previous 
one (12] . 

It is proved that only the "Kutta condition" solution ful- 
fills the boundary conditions of equal displacement and 
pressure on both sides of the shear layer. Within the 
limits of this model, it is proved, that only forced in- 
stability waves exist. It is shown, that no ambiguity is 
left in the excitation mechanism. 

For the parabolic excitation field, also analytic solutions 
for u^ and u 2 at both sides of the shear layer are calcu- 
lated. For x > o this solution consists only of a damped 
and an amplified instability wave. 

Numerical calculations of u 2 (velocity fluctuations in the 
region of existing mean flow) are carried out for y 4“ o. 

The u 2 computations are of crucial importance for the 
comparison with experimental data. The calculation utilizes 
a source distribution approach and is based on the analyti- 
cally given v 2 velocity distribution. Also simple 
asymptotic formulae for the field outside the interaction 
region at the plate edge are given. Thus, the region of 
numerically computed velocities is completely imbedded 
into asymptotic equations and fits also the analytic 
solution at y = -o. In addition, the computations have 
been carried out by two different persons using different 
mathematical procedures, computer languages and computers. 
Therefore, the numerical computations can be considered as 
a reliable bridge between the analytical solutions for 
u and v, derived from separate calculations. This is also 
part of the proof 2) that the forced instability waves do 
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fulfill all boundary conditions and additional arbitrary 
instability waves do not exist. 

5) The excitation by different types of sources in various 

locations is discussed. With the exception of the (singular) 
cases if either 

(a) the excitation field has a stagnation point at the 
plate edge or 

(b) the excitation is directly located at the plate edge 
itself# 

the exteriour excitation creates a parabolic pressure near 
field at the plate edge which dominates the control of the 
shear layer. The possible existence of a stagnation point 
(case (a)) is demonstrated with a dipole excitation at a 
certain orientation with respect to the plate. The effect 
of a monopole source on the shear layer is investigated in 
detail (see Fig. 21 ) . 



Fig. 21 Excitation by a monopole source in different 
locations 

A monopole in position (a) nearly creates a stagnation 
point at the plate edge. Therefore, this type of excitation 
is very inefficient. The excitation by a feedback from down- 
stream turbulence would be also of this type. A source in 
position (c) would produce a strong parabolic excitation 
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field and would be by orders of magnitude more efficient 
than an excitation in position (a) . This is the reason why 
free shear layers are so sensitive to sound waves from up- 
stream and so insensitive to the strong near field pressure 
created by the turbulence. 

The excitation directly at the plate edge (b) has been 
used in some recent experiments. It was not clear, however, 
what the relevant reference quantities were. An expansion 
for small distances provides an elegant and simple result: 

u(x-|y[ ) 

n 5 0 

<183) |u,| = %- • e for — > 1 

* 2U U 

' o o 

For this case also a complete solution for the shear layer 
motion at all x can be given, as well as an analytic 
solution for the whole induced field at y *j* o, in terms of 
exponential integrals of complex arguments (this function 
is tabulated in [31)). The result in eq. (183) refers to 
the fact, that 1/2 of the volume flux Q of the source pene- 
trates through the line y = o. In a real experimental 
situation, this ratio might be higher, maybe closer to 1, 
depending on the individual configuration of the excitation 
device at the lip. 

6) The major part of this paper, however, deals with the 
excitation by a parabolic field at the plate edge, which 
is created by any source far away from the edge. It is 
shown, that in this "parabolic" case the interaction takes 
place within less than one wavelength (of the instability 
waves) downstream of the plate edge. Downstream of this 
interaction region, the v-velocity field (at y = ±o) splits 
into one contribution of a pure amplified instability wave 
and another contribution from the excitation alone. 

7) The theory is also extended to the "two stream case", where 
we have two different mean velocities and densities above 
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(0^ , p ^ ) and below (U 2 f P 2 ^ s hear layer. Although the 
details of the calculations become somewhat awkward, the 
mathematical results are still very similar to the preced- 
ing ones. The quantity in the excitation field that is 
most easily dealt with is again the pressure, and it 
becomes clear, that pressure gradients (3p/8y) rather than 
velocities drive the shear layer. Again, the far field 
downstream of the interaction region splits into an 
amplified instability wave and a forcing velocity distri- 
bution which is continuous through the shear layer. Since 
the densities on both sides of the shear layer are 
different, nature prefers the average value of both densi- 
ties as a coefficient in front of this forcing velocity. 
The consequences of the introduction of the new parameters 
o = U-|/U 2 and 0 = p^/p 2 are the following: 

(a) the downstream amplification rate of the instability 
waves is reduced for reduced difference of the veloci- 
ties, as expected; 

(b) if the low speed region has a lower density than the 
high speed region, the amplification rate is also 
reduced. 

The present theory predicts potential fluctuations outside 
an infinitesimally thin shear layer. In an experiment, 
however, we have always a shear layer of finite thickness. 
Nevertheless, at low Strouhal numbers S„, both theories 

U 

provide the same potential field outside the shear layer 
(if the instability wave is fully established downstream 
of the plate) . By matching both (identical) potential 
fields, we can predict also the motion inside the shear 
layer, as calculated by Michalke [6 ] . 

Obviously, the influence of finite shear layer thickness 
leads also to deviations from the present model. If the 
data acquisition in an experiment relies on measurements 
in the potential field outside the shear layer, we have 
to consider in how far the wave number in the real case 




deviates from our idealized case. The wave number is an 
essential quantity, because it determines the decay rate 
of the induced field perpendicular to the shear layer 
(y-direction) . Prom Michalke's [6] and Freymuth's [7] 
results it is clear that the wave number becomes higher 
than in our model for increasing S^. This leads also to 
an increased decay rate in the y direction. 

10) Another typical deviation in experimental setups is the 

"overshoot" (mean flow gradient) near the shear layer. It 
is not too difficult to model this. The analysis yields 
an increased wave number, which again enhances the decay 
rate of the induced velocities in the y-direction. Inter- 
estingly enough, the downstream amplification rate of the 
instability waves remains almost unchanged. 


11) In a typical experimental situation the two streams at 

both sides of the shear layer are limited in their width. 
Thus, in reality, we have a shear layer in a channel 
rather than a free shear layer in an unlimited flow. 

There are several things to be considered: 

(a) the general solution of the shear layer motion is 
still the same, 

(b) the excitation pressure field, however, is changed. 
There is no difficulty, though, to calculate it with 
conformal mapping. The excitation pressure field is 
still parabolic in the vicinity of the splitter plate 
edge, but it falls off exponentially beyond a certain 
distance in the downstream direction. 

(c) With this new excitation pressure field the shear 
layer motion can be calculated, at least downstream 
of the interaction region. The analytical solution 
contains r-functions of complex argument. Since this 
solution is quite awkward to handle, it is expanded 
for large md/0 o , where d is the channel halfwidth. 
This means that the wavelength of the instability 
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waves should be relatively small compared to the 
channel half-width. We end up with a reduction of the 
instability wave magnitude, represented by the co- 
efficient (1-nwd/160 o ) . In comparison to this quite 
significant influence on the magnitude, it turns out, that 


(d) the direct reflection Influence of the channel walls on 
the velocity distribution is quite insignificant (close 
to the shear layer) as far as the instability wave 
downstream of the plate edge is concerned. 


12) Perhaps the most significant finding of the present paper 
is a reference quantity for the excitation of a shear layer 
by a parabolic pressure field, which is a fairly general 
ease, as discussed in paragraph 5). It is the 


reference quantity Ap 12 //1 


The significance of Ap^ 2 anc * 1 can be seen in Fig. 22 . 
ap 12 //l is related to the excitation field on ly. It will 
not change, e.g., if the mean flow U Q is switched off. For 
instance, the u 2 -velocity of the instability waves close to 
the shear layer depends on ftp^ 2 //l in the following way: 


(184) 1 u r 


!&P 

7 ? 


1 2 1 . 1 


p > / wU" 
o 


/tT 

4-V2 


valid for wx/q > 2 


89 







\ 

< 




ORIGINAL P*®*** 
OF POOR QUALITY 


Ap, 2 ^.excitation Ap L _^ -v^exciran 

.. , y J / J--I 

— 7 Uo 1-H ^ - .- » u »y ^ 


excitation 


Excitation in free Excitation in a 

conditions channel 

Fig. 22 Configurations for shear layer excitation . 


Equation (184) is well confirmed by recent, yet unpublished, 
measurements of the author. 


For an experiment in a channel (see Fig. 22), the situation 
changes slightly. If the two microphones used to measure 
Ap at both sides of the splitter plate are moved further 
upstream (to improve the accuracy of the measurements) , we 
have 


(185) 


Ap 12 _ ajd . 2 /v 

ft ^ , iu 

In — * — » — 

.i- 



The significance of L and d can be seen in Fig. 22. The 
magnitude of the instability waves is then the value pro- 
vided by eq. (184), but multiplied by the coefficient 
(1-7rwd/16 0 Q ) . The latter coefficient (see paragraph 11 ) 
reflects the limited size of the parabolic pressure distri- 
bution at the plate edge in the channel. 
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10. Appendices 

10.1. Ap pendix A 

Glossary of symbols 

a Q sound speed 

a 1 ia^ constants in Urn numerical representation of the 

error function erf of a real argument 

b^ , bg lower and upper boundary of a vorticity distribu- 

tion on the x-axis (see Fig, 7) 

phase velocity of an instability wave 

rel ative phase velocity of an instability wave 
c ph rol a C ph/U 2 

d half width of a two-dimensional channel (see 

Fig, 17) 

z 

2 ! - 1 ’ ^ 

erf(z) error function, defined as erf(z) * < — e dt 

” /7t i 

o 

orfc(z) complementary error function erfc(z) <* l-orf(z) 

f frequency of the excitation f * 

g^,g 2 lower and upper boundary of a source distribution 

on tho x-axis (see Fig. 8) 

h displacement of the shear layer duo to the insta- 

bility waves and the excitation field 

1 imaginary unit V » *'-T 

1 distance between the edge of the semi-infinite 

plate and the pressure pickup devices (e.g. , 
microphones) at both sides of the plate 


c ph rel 
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m, n summation variables (integer) 

p pressure, usually the perturbation pressure , i.e., 

the deviation from the ambient (atmospheric) 
pressure. All perturbation quantities are pro- 
portional to e 

p 1 fluctuating pressure above the shear layer at 

y > o 


P 2 fluctuating pressure below the shear layer at 

y < o 

Ap, AP 12 pressure difference between both sides of the 

semi-infinite plate 

p abbreviation of a quantity related to a pressure 

gradient in y-direction, i.e., 

. _ 2i . 3p 2f 

” 3y 


p if' p 2 f 


contribution of the exteriour forcing field to 
the total fluctuating pressure. P-]f = P 2f 


Pis' p 2s 


q, q ( 5) 

r 



contribution of the shear layer field to the 
total fluctuating pressure. p^ s is valid for 
y > o and P 2 S is valid for y<o. In addition we 
have Pl ■« p 1f + p 1s ; p 2 - P 2f * P 2s . 

source distribution on the x-axis> dimension like 
a velocity. 

radius or distance from the origin (in most cases 
the splitter plate edge) to the considered field 
point r =/ x 3 + y s 

distance between the location of a point source 
and the splitter plate edge r D = /x Q J + y Q 2 
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S3 


integration variable (eg. (38)) 


t time 

t integration variable, e.g., in oqa. (171) anti 

(173) 


u fluctuating velocity component in the x-dlreetlon s 

all fluctuating quantities are proportional to 
Q -.Ut 

u., fluctuating u velocity component above the shear 

layer at y > o 

U|^.| u- velocity induced by a v velocity distribution 

of an amplified Helmholtz instability wave 
(exponent!,) in the half-plane above the shear 
layer (y > o) , The v distribution of the instabi- 
lity wave exists only for x > a, 

similar to u,, jV, but referring to the attenuated 
Helmholtz instability wave (exponent \ 2 ) 

fluctuating u- velocity in a channel with a 
splitter plate, asymptotic value far upstream of 
the trailing edge of the splitter plate 

fluctuating u velocity component below the shear 
layer at y < o 

similar to u^., » but valid below the shear layer 
(y < o) 

v ' 2 ,\ 2 similar to u, V2 , but valid below the shear layer 

(y < o) 

u 2l) u-velocity , induced by v^ at y < o. v ^ is the 

difference between the forced instability wave v a 
and the asymptote which behaves like the Helmholtz 


1 -V2 


jsasasKsafERB?^ 
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AU 2D 

U 2DI 

U 2DR 

U 2I 

U 2R 


instability wave v 2H ; v 2D *> v 2 - v 2n . v 2p becomes 
zero both at x «* - » and x ** + <», 

u-volocity, induced by just one trapezoidal source 
element representing the local v 2D distribution 
between two neighboring points £ a and 

imaginary part of u 2D 

real part of u 2D 

JLmag inary part of u 2 

real part of u 2 

u^-velocity component induced by a point source 
located at y = o and x = o. 


u c L center line velocity (at the splitter plate) in the 

two-dimensional channel 

u velocities (with 
above) in the following 


in the two stream case P has to be replaced by 
P 2 and U Q by U 2 

u ind velocity ind uced by the mirror image of the shear 

flow representing the wall in a channel flow 

v fluctuating velocity component in the y-direction, 

all fluctuating quantities are proportional to 
— iuifc 
e 

v 1 v-velocity component above the shear layer at y > o 
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u dimensionless form of the 

various indices, as given 
form 

u • p el 
u - 

6p,j 2 
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If 


'Is 


'ID 


II 


1 R 


1DI 


1 DR 


v-velocity component at y > o, generated by the 
extoriour forcing field 

v-veloaity component at y > o, generated by the 
shear layer; v 1 * v 1f + v 1g 

difference between the forced instability wave v^ 
and its asymptote v^ H , which behaves like an 
amplified Helmholtz instability wave, v^^v^-v^j 
(valid for y > o) 

imaginary part of v. (valid for y > o) 
real part of v^ (valid for y > o) 
imaginary part of v 1D (valid for y > o) 
real part of v 1Q (valid for y > o) 


III 


asymptotic behaviour of the v^ component of the 
forced instability wave, assumes the form of a 
Helmholtz instability wave, amplified in the down- 
stream direction (valid for y > o) 


2D 


2DI 


2DR 


2H 


v-velocity component below the shear layer at 
y < o 

difference between the forced instability wave v 2 
and its asymptote v 2H , which behaves like an 
amplified Helmholtz instability wave. V2D =V 2“ V 2H 
(valid for y < o) 

imaginary part of v 2D (valid for y < o) 

real part of v 2D (valid for y < o) 

asymptotic behaviour of the forced instability 
wave, assumes the form of a Helmholtz instability 
wave, amplified in the downstream direction (valid 
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dimensionless form of the different v-velocities 
as given above, i.e., 

v • pi's Ml 

v = t— — ; in the two stream case p has to 

P 12 

be replaced by p 2 and 0 Q by 0 2 


v-velocities at the lower and upper end of a 
trapezoidal source distribution. The local 
strength of a plane source distribution is two 
times the local v-velocity just above the distri- 
bution 

horizontal coordinate 

nondimensional horizontal coordinate x=xw/0 or 
o 

x-xw/Ug for the two stream case 


vertical coordinate 


nondimensional vertical coordinate 
y=yio/u o 

complex variable, usually defined as z=x+iy=re il ^ 1 
complex conjugate of zj z = x-iy = re - ^ 


constant, used on different occasions: 

in eq. (31) it defines the strength of a pressure 

field 

in eqs. (67) and (140) it defines the strength of 
a potential 

in eq. (155) it is a constant with the dimension 
of a velocity, A = u^ 

in the appendix D, eq. (D3) , A is the real part 
of a complex variable in conjunction with a 
computer program which uses the same symbol 


* «« 


ABS(z> 




, < 
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modulus of a complex quantity , e.g., 
ABS(z) * /x 3 + y J 


B 


B 


B (x,y) 


C 1 ' C 2 


D 


E 1<z) 




constant, used on different occasions: 
in eq. (67) it defines the strength of a potential; 
in eq. (160) it defines the strength of a pressure 
field 

in the appendix D, eq. (D3) , B is the imaginary 
part of a complex variable in conjunction with a 
computer program which uses the same symbol 

Beta function, Euler's integral of the first kind 
(see also [31 , 39)), defined as 

1 

B (x,y) - i 
o 

constants (can be complex) with the dimension of 
a velocity, defining the magnitude of the ampli- 
fied (C^ ) and decaying (C 2 ) Helmholtz instability 
wave 

Dipole strength 

Exponential integral, defined as (see also [31}) 


00 



z 


lower boundaries of the integrals in eq. (23) ; 
their existence is equivalent to the existence of 
additional Helmholtz instability waves 

lower boundaries of the integrals in eq. (66) 
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H (x) , Uoovisido unit stop function llj x j B o for x < a 

and It * 1 for x > o 

I imaginary part of a complex quantity as used in 

the appendix D, z *> R + il, in conjunction with 
a computer program using the same symbols 


imi I 




M 

o 


integral, defined by eg, (169) 

imaginary part of a complex quantity in brackets 

distance between splitter plate edge and the 
location of the pressure pickup devices (e.g., 
microphones) at both sides of the splitter plate 
in a channel (see Figs. 19, 22) 

Mach numor, e.g., M => 0 Q /a o 

«• 2 

e.g., in Q(z ); of the order of z 


P 


Q, Q, 


R 

R 

R 

Re{ } 
RS 


Fourier transform of the pressure 

points above and below the shear layer at which 
induced velocities are considered 

source strength Q « Q e” 1 '*'* 1 ; Q,Q o have the dimen- 
sion of a velocity times a length 

radius of a nozzle (used only in the introduction) 

cylinder radius (eq. (31)) 

real part of a complex quantity as used in the 
appondlx D, z - R + iX, in conjunction with a 
computer program using the same symbols 

real part of a complex quantity in brackets 

right hand side of eq. (119) 
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source strength of a pressure source 

Strouhal bomber calculated with the momentum 
thickness of the shear layer and defined as 
s tt ■=> f * n/u 0 

abbreviations for complex quantities 
i (r. »x ) _ Urtx) 

T « 1 t _.^ : 0 J ,,, a 1 „ 

* o J o 

mean flow velocity, it has only a component In 
the .'{•‘direction 

mean flow velocity for y < a 

mean flow velocity (in the two stream case) for 
V > o 

moon flow velocity tin the two stream ease) for 
y * o 


angle of the radius vector in the a ••piano} 

•f ft X + Jv e r O’*’' 
a a a 


complex wave number as defined usually in the 
stability theory. The relation to our nomenclature 
is the followings l,v *> \ 

angle of the radius vector In the s, -piano, 

hi 0 * 

z b & *b * e V 

ratio of the densities at both sides of the shear 
layer.; 1? « ey/e^ 

strength of a vortioity distribution, located on 
the x-axis, It has the dimension of a velocity 




delta function at the location x * x 
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quantity defining the slop'* of a shear layer 
C * 1/29 


* 


Fourier transform variable, corresponding to y 


angle of the radius vector in the z-plane; 

. l*S 

z * x + iy « re # 


complex wave number of the instability waves in 
a shear layer. The relation to the conventional 
wave number a in the stability theory is x * ia 


' 1,2 


Simplex wave number, defined 

U) for a flow below the shear layer with U*U 0 
and zero flow above the shear layer as 


' 1,2 


U 


(1+1); the positive sign corres- 


ponds to X ^ and the negative , v n to X 2 

for a flow with different velocities (0^ , 
U 2 > at both sides of the shear layer and 

, 0 2 ) as 

l(1+o3) *■ ( 1 — a ) * /lj 


with different densities (p 1 as 


- HL . 1 

l 1,2 r . 1+a 2 3 

U 2 


< iii) for a flow below the shear layer with uni- 
form gradient du/dy and the value 0s0 o at 
the shear layer (y--o) and with zero flow 
above the shear layer 


- JL 

‘1,2" 0" 


in* 2 


1,2 


is the corresponding dimensionless wave number 
with x. „ = x. _ • u o or x. „ * U 2 for the two 

\ t £. \ f C I 9 *• “T* 

stream case 
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Fourier transform variable, corresponding to x 

integration variable for source and vorticity dis* 
tribution on the x axis; 5 has the same direction 
as x 

lower and upper boundaries of a trapezoidal source 
distribution element 

lower boundary for the integral over a continuous 
source distribution with q ^ * 1//lf 

upper boundary for the integral over a parabolic 
source distribution q^ with q^ « /? 

stepwldth (length) of a trapezoidal source element 

density of the fluid 

density of the fluid above the shear layer (y > o) 

density of the fluid below the shear layer (y < o) 

velocity ratio of the two stream case with 
a* U 1 /U 2 , 0 being the mean flow velocity above 
the shear layer (y > o) and U 2 being the mean flow 
velocity below the shear layer (y < o) 

complex constant in eg. (93) 

amplified and decaying Helmholtz instability waves 

radian frequency of the exteriour excitation 
u = 2 ii f 

Gamma function) Euler's integral of the second 

kind, defined as 

f 00 

r (z) = { G_t * fcZ ' 1dt 
o 
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0 momentum thickness of the shear layer, defined as 



<|> potential function for tho fluctuating velocities 


<f> stream function for the fluctuating velocities 


[£) complex potential function, defined as $ * $ + 

(/* Laplace operator in two dimensions, defined as 

y l p * + -jy? for p, as an example 


10.2, Appendix B 

Discussion of the pressure sources 

The source term in eg. (11)', on the right hand side 


( 11 ) 


v J p 


•2p 


DU 

Dy 


_3v 

8x 


<* S 


has apparently led to some confusion, in particular among those 
people, who have heard about aeroacoustic theory and know for 
sure, that there should be "some quadrupoles" [40] . We will not 
query Lighthill's theory, but the present approach is something 
different. We will have to deal with pressure sources? these are 
not sources of matter, which would be at variance with the 
continuity equation in our model (3). 


To provide a better feeling for our pressure sources we will de- 
termine the pressure source structure of an instability wave in 
a free shear layer. This refers to the right hand side of eq, 
(11). Using a Fourier transform technique, we integrate eq. (11) 
to obtain the pressure distribution in the whole field, which 
could be obtained clearly also by other methods. Finally, we 
check our solution with the Euler equations. 
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Wo restrict ourselves to the "one stream" case where we have a 
fluid at rest above the shear layer and a flow of velocity V Q 
below the shear layer. Wo consider a normal spatially amplified 
instability wave where all fluctuations like p, v, u and the 
elevation h are proportional to j n addition, we restrict 

ourselves to low Strouhal numbers (calculated with the shear 
layer thickness) whore wo have \ * (i+1) *a>/U * It turns out to 
bo most convenient to replace the v-velocity in eq. (11) by the 
displacement h. The displacement (as the pressure) has the same 
value at both sides of the shear layer and can be assumed to be 
constant throughout the shear layer at low Strouhal numbers. The 
relation between h and v is 


( 7 ) 


111 + 0 
at u (y) 


.ix 


inserting this into eq. (11), wo have 
3 dU 


(B1 ) 


S ° - 


2 0 co a 

u 

o 


All 


h 


dy 


(k,y,t) 


(1 -1(1-2 -M)) 


by writing h or h| x ^ wo mean h e ^r-iuit throughout the cal- 
culations. It is assumed, that h is constant in the shear layer) 
it will vary, however, for greater distances |y| from the shear 
layer, as wo will see later. 


A reasonable model for the shear layer velocity distribution is 
a hyperbolic tangent profile [6, 7] s 


(B2) U 


u 


(y) 


~ (1-tanh r,y) 


The quantity 5 in this equation is related to the momentum thick- 
ness 0 


•f M 


(B3) 0 » 


_lill (1 lid)dy 

u_ u. 


Ill 






\ 
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* % 




in the following way 
(B 4 ) c * 1/20 
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Fig. 23 Velocity distribution in the shear layer 

Using eq. (D2) for the velocity distribution* we obtain the real 
and imaginary part of the pressure source distribution S 

1 

cosh 1 (y/20) 1 

tanh (y/20) 
cosh J (y/20) 


The distributions in y-direction of these functions in the shear 
layer are displayed in Fig. 23. 

As we see from Fig. 23, the real part of S has the near field 
structure of a S-f unction’, it exists only in a limited region 
which can shrink to zero if Q -► o. The actual value of the 
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& function is determined by the Burface under Uho 1 /cosh’ -curve 


in Fig. 23. 


(D6) | oo^rrnr * 40 ■ 


Therefore, we can replace s Ro by a 5 function representation for 
0 * o 


(07) S Rq - 2pu>»h 5 (y) , 

Tlve imaginary part has no contribution of a 6 function typo, and 
an integral like (n<5) would give zero. Contributions of a <$'- 
type (dipole) are pirosent, but this influence becomes negligible 
with o ► o. Therefore, the contribution of the imaginary part of 
S becomes zero. We are loft with the nonhomogeneous partial 
differential equation 

(08) V* p * 2pu , h 5 (y j , 

This equation will be solved by Fourier transformation. In fact, 
using h and p proportional to Q^x-imt wc j iave already taken a 
Fourier transformation in time (harmonic motion) and in the down- 
stream direction x. To obtain a notation compatible with the 
usual form of Fourier transformations wo replace \ by iv , 
whore v is the Fourier transform variable corresponding to the 
physical variable x. Wc have, therefore, Instead of eq. ( B8 ) : 

(B9) -v 5 p + * 2pw’h $ <y ). 

We take the Fourier transform in y-direction‘, with 

•fo> '*•<« 

( “ 101 p („> ■ j ol " V • e|y)‘ ,yl p (y> ■ A/ °" lnY p (,i) a ’- 

-»i “TO 

(BID — v 3 P-n J P ■ 2p« a h 
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where h is the displacement of the shear layer at y * o. The 
transformation back into the physical plane can be found in 
tables of Fourier transforms, like C 39) 

<B13) p - - 


For h we have 
(B14) h 


Xx 


= h • e nA with iv = X = — (i+1). 
(y-o) °o 


This can be inserted into (B13) and we find 


(B15) 



The last two terms can be interpreted as the spatial distribution 
h(x y) t * le ^ is Pl acement 

(B16> h (x v» - • e <VU 0 )-(i + 1)-(x + i|y|) 

lx# yf o 


and we find the relation 

(B17) 1 ** (x»y) 


(Hi) 


pwU o ’ h (x,y> 


We can check whether or not this relation is correct by using the 
Euler equations and the relation between v and h, i.e . , eq, (7). 
For y > o, in the fluid at rest (index "1"), we have 


V 1 ' -iu)h (x,y) * 
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Hi 


4 


\ 

4 * 


With the first Euler equation, we obtain 
(819) ^ ■ f'Vyl* 
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If we take the y-derivative of eg. <B15) we find, that this 
relation is fulfilled. 


With that result, it is shown, that it is possible to analyze 
the pressure source distribution within a shear layer and inte- 
grate it (via Fourier transform techniques) to obtain the 
pressure distribution in the whole flow field. Thus, there is 
nothing mysterious left in our pressure source approach. 

We have to mention what the difference between a pressure source 
(S) and a volume source (Q) is. A pressure source is a singulari- 
ty in a field which is described outside the singularity by the 
Laplacean equation 7 J p = o. A volume source is a singularity in 
an analogeous field with V 2 <j> = o. In a fluid at rest, the re- 
lation between both fields is simple. From the first Eulerian 
equation we have iwp-|^ - i.e., both fields are proportional. 

This is the reason why we have used potential theory and veloci- 
ties like the excitation velocity v^ in several sections of 
this paper instead of using pressures and pressure gradients. It 
would be more appropriate, however (and less obvious), to use 
pressure sources throughout the whole paper. The basic differ- 
ences occur in the shear layer. Clearly it makes no sense to 
define a velocity potential inside the shear layer, whereas the 
pressure is still well defined there. In addition, the existence 
of pressure sources in the shear layer does not violate the con- 
tinuity equation. 

Some readers may be dissatisfied with the above methodic approach- 
Indeed, we have used a right hand side of eq. (11) which is known 
to fulfill all our equations. Instead of dealing with a formally 
nonhomogeneous equation we can also completely eliminate v and 
Write eq. (11) only in terms of the pressure. This can be done 
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by using the second Euler equation. We end up with 


(B20) (~iu> + D = 2 


a*p 

3x3y 


outside the shear layer we have again f'p « o, If we insert an 
instability wave ansatz p « e^ 0151 we find 



This is an equivalent to the Rayleigh equation for the instabili- 
ty of inviscid parallel flows, written in terms of the pressure. 
Outside the shear layer we have p w e a ^ , as usual in the po- 
tential field of a wavy motion. The complete distribution of p 
can be found using a numerical procedure similar to the one used 
in stability theory [6] . This procedure would provide results 
valid for shear layers of finite thickness and would produce the 
whole pressure field for an instability wave, extended from 
x = ■» to x = +<■>. In fact, this latter approach provides a 
genuine and complete solution of the pressure field without 
knowing anything in advance, whereas the approach with a non-^ 
homogeneous pressure equation requires previous knowledge of the 
fluctuations in the shear layer. 


10.3. Appendix C 


The absence of "free" instability waves 

In the main part of this paper, forced instability waves are 
calculated, which is equivalent to the assumption that the 
lower boundaries of the general solution of the instability waves 


B dx 


, 23 ) v t u-t e X 1 x [ e'MV^dx + § e X 2 x J e" X 2 x v 1f c 

O p P p 

G 1 G 2 

should be - o. This would also correspond to the non- 
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existence of additional "free" instability waves of the Helmholtz 
type. In a previous publication of the author 112), it had been 
mentioned, that additional instability waves of the Helmholtz 
type do not fulfill the boundary conditions, but no detailed 
proof had been given there, it is quite obvious, that a solution 
valid for a sheaf layer extended from x * -» to x * +«* cannot be 
valid for a semi- infinite shear layer with a rigid boundary, such 
as a semi-infinite plate. The reader who is convinced by this 
simple and conclusive argument does not need to read the remain- 
der of this section C of the appendix, 

However, with all the present confusion about the Kutta condition 
issue the author feels that he cannot avoid providing a mathema- 
tical proof for the fact, that the Helmholtz solutions do not 
fulfill, in any conceivable combination, the boundary conditions. 
The absence of Independent Helmholtz solutions leads to the con- 
clusion, that only forced Instability waves exist and that only 
the Kutta condition is possible . Of course, this latter state- 
ment on the Kutta condition is only valid for an infinitesimally 
thin shear layer, as assumed in our model. In reality, this would 
correspond to the case of low Strouhal number S Q . Any "rectifica- 
tion" effect (i.e., a steady deflection of the jet caused by the 
Unsteady motion) as suggested by Orszag & Crow [9] could not be 
found in the experiments [41,42]. Such an effect would be also 
at variance with the equations of motion for a linearized 
problem. The linearization is well confirmed by the experiments 
(41, 42). 

In the following calculations, we will restrict ourselves to the 
one-stream model of section 2, To prove that free instability 
waves (Helmholtz solutions) are not solutions of the present 
problem, we have to keep in mind that these solutions are 
truncated for x < o, because v has to be zero at the semi-infi- 
nite plate. The v distribution, therefore, assumes the following 
form: 

(Cl) v 1 = H (x) (C^* 1 * + C 2 e A 2 x ) 
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In this equation is the v- component of the fluctuating veloci- 
ty just above the shear layer at y = +oj C 1 and C 2 are constants 
which may be complex and Hj x j is the Haeviside unit step function 
which is zero for x < c and unity for x > o. From equation (Cl) 
we can calculate v 2 , the v-velocity just below the shear layer, 
using equation (10), which represents the condition of equal 
displacement at both sides of the shear layer: 

noi v 2 - .,.1^ 


we obtain 


(C2) v r 


■ iH (x) [v* 1 * - c 2 el2X ] * ‘TT 


[VS 


The delta function & is created by taking the derivative of 
the step function H | x j* The physical significance of the 6 ^ x j 
function is that we need a source to create a finite displace- 
ment of the shear layer at x » +o in the flow (at y < o and with 


mean flow U Q ) . 


After having fulfilled the condition of equal displacement at both 
sides of the shear layer, we have to check whether or not the 
pressure is the same at both sides. Instead of taking the 
pressure, we can take as well 3p/3x, the pressure gradient in 
the downstream direction. 3p/3x is related to the u-velocity 
component via the first Euler equation. We end up with eq. (6) 


u, + i— 

2 oj 


i 3(prp 2 > 


The right hand side should be zero for x > o. We have to calcu- 
late first the induced u-velocity for the decaying Helmholtz 
wave (exponent X 2 x) . In addition, we can omit the Haeviside step 
functions, if we change the boundaries of the integrals 


2tt J z- 


X 2 5 C 2 r 

“ d? + 57 j 


e X 2? 
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It is understood, that the whole induced velocity u^ will then 
consist of 


<C4) u, - U,^ + u U2 


The integrals in eg, (C3) are of the structure of exponential 
integrals, defined as (see [31]): 


(C5) E 


-t 


1 (z) 


e~ 

t 


at . 


After appropriate substitution of variables we obtain 

C, 


,c6 > “n, * - if P 2 " • W> * e>25 • E i'V>] • 


This is the induced u-velocity for y > 0 and for an instability 
wave assaying in the downstream direction. The situation is 
slightly more complex for an amplified instability wave (terms 
containing By a straight forward integration we would 

enface divergence problems with the upper boundary of the inte- 
grals. This problem can be circumvented if we subtract from the 
complete instability wave solution (extended from x = to 
X = +») a contribution from x = -« to x = o to obtain v^ = o for 
x < o. The relation between u and v of the complete wave is 
(see eq. (68)) - -iy^ , Otherwise, the integration is very 

similar as that leading to eq. (C6). We obtain after some inter- 
mediate calculations 


(C7 ) u = -iC^ 1 * - [^ A 1 Z • E 1 (i 1 z)te A 1 z • E, (^zjJ . 


The complete solution for u 1 is then (y > o) 


(C8) u. = -iC.e A 1 


ilX 


^ | I.— 0mm 

- [e A 1 z ■ E 1 U lZ )+e A 1 Z • E 1 < A 1 5 ) 


- ^ P 2Z ‘ V A 2 z)+ eX2Z • E 1 (A 2 Z) ] 

for y++o . 
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It is understood, that y becomes +o to obtain the valuo of u^ 
just above the shear layer. This transition to y • +o can be 
postponed, though, to a later stage of the calculations. 


The procedure to obtain u 2 is very similar. However, for y < o 
the sign of u 2 changes, in addition, we have a contribution from 
a point source, represented by a 5 function appearing in u,. For 
negative y we have as the contribution from the source term 


<C9) 


u 26 w tt ] <x-5) * + y 2 d5 

— w 


For y *> -o we obtain 
iU 

' c, °> “ 2 , * - iS <V° 2 > ' 


Obviously, the induced u-field of a source located on the x-axis 
should be proportional to 1/x. The complete solution for u 2 is 
then from eq, (C2) 


(Cl 1 ) 


u 2 = -C 1 e A 1 X +i ~ Je X 1 2 - E^^zJ+e^ 2 • E 1 < X 1 z ) j[ - 

- sr E* 2 ' • e 1 <V M * 2 ’' E A ! >] - «V c a> 


for y-t— o, 

u,j and u 2 have to be inserted into equation (6) to find out 
whether or not the pressure equilibrium at both sides of the 
shear layer is fulfilled for x > o. To do this, we have to take 
the derivative of the exponential integral 


(Cl?) 


dE 1 ( z ) 
dz 



Using this relation and, after taking the derivatives, letting 
y o, we obtain from equation (6) after some algebra 


(Cl 3) 



IT tDX 


(<c 1 -c 2 l-i<c 1 tc 2 l>- ( Cl .c 2 ) 


j_ d( PrP2 1 

pw dx 
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The right hand side should be zero for x > o. We see, that the 
right hand side can be zero for nil positive X only if c^“C2*o. 


This is the proof for the absence of "free" (unforced) instabili- 
ty waves- euu it is, since the -orced instability waves fulfill 
the Kutta condition at the plate edge? the proof for the existence 
of the Kutta condition there. Here- the Kutta condition is under- 
stood in that way, that the flow does not suffer a local jump ill 
its slope when it passes the trailing edge of the plate. 

10,4, Appandlx D 


v-computatlons at y * <o, tables 


The purpose of this section is to provide numerical data for v,. 
and v 2 - This means, that the error function of a complex argument 
has to be computed. The present calculations are based on ex- 
pansions given in (31], One of the building blocks of our cal- 
cuxation is to provide an expansion for the function e erf z 
(Where z is an arbitrary complex variable) and to split it sub- 
sequently into real and imaginary constituents. In the preceding 
calculations wc used this expansion; 


( D 1 ) o’ 3 erf z 


2 

77 


£ 

n=o 


1-3. . ."'(Sn+1) 


2n+1 


It provides, e*g., the near field behaviour near the plate end 
in a vary elegant way. It was also used to produce the numerical 
values for v 1 between x « o and x - 2 computed in (12). Unfortu- 
nately, numerical convergence problems occur beyond x a 3 so that 
this serias cannot be used up to those values of x where the 
asymptotic expansion for large x (say, x > fij 


(D2 ) 


5»J 

e 1 erf z 
(z •» «) 



tfJ 

1+ 2 (~1)‘ n 
m-1 


1-3... (2m-i ) | 
(2z s ) m > 


is valid. The series expansion in eq» (D2) is semi-konvergent. 
Therefore, it does not make sense to include, say, more than 
6 terms of the series in our case. 
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For our numerical computations wo had to use a mote complicated 
expansion providing numerical convergence for low and middle 
values of the argument so that also the lower validity range of 
eq. (D2 ) (for large z) can be Covered. This is based on eg. 
17,1.20) in (31). We introduce z * R + iX and 

<D3> A + iB * e z> erf % « e (R+iI) a erf (R+il) . 


After splitting eq. (7.1 .29)ln (31) into real (A) and imaginary 
parts (B) we obtain after some algebras 

_ 

<D4) A * e R _I ■ coe2RI • erf R + |^- (cos 2RI-1) + 


2e 


-i ? 


n* 

"T" 


E -TTXrrr (2RCOS2RX ~ 2R cdsh nl) 
n *= 1 n +4R 


(D5) 


B * e 1 * 3 " 12 . sir,2Rl • erf U + • sin 2RI + 


-I* 


+ 2S-1- s (2Rsin 2RI + n sinh nX> 

n=1 


For the numerical representation of the error function of a real 
argument we use an approximation suggested by Hastings (31 , 
eq. (7, 1.28)1 ( 

(d 6 ) erf R s i - « — 1 — — 

(1+a 1 R+a 2 R , +...+a (3 R 6 ] lc 

with the constants a : 
n 

a 1 = 0.0705230784 a 2 * 0.0422820123 

a 3 = 0.0092705272 a 4 = 0.0001520143 

a 5 « 0.0002765C72 a g * 0.0000430638 . 

The error of this approximation should be smaller than S^IO - * 7 
over the whole range. 
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With equations (D2) and (D4), (.95) we can now write the computer 
subroutines for the building block "e z erf z" which is the 
essential part for the v 1 and v 2 computations. When we will show 
some computer programs here, it is done without showing all de- 
tails of the considerations leading to these programs, The pro- 
grams shown here are written in BASIC for the Bharp PC 1211 
Microcomputer* > (see Pig. 24). These programs are only provided 
in order to keep the mathematics verifyable. In addition, it is 
shown, that the present numerical solutions are not dependent on 
the access to big, fast and expensive computers. 

In the subroutines, which are given below, R and I are, as in 
eq, (D3), the real and imaginary part of the variable z, and 
A & B have also the same meaning as in eq. (D3) , For equation 
(D4), (D5 ) we have the subroutine "SUBR 1" and for equation (D2) 
we have "SUBR 2". With these subroutines we have to calculate v 1 
and v 2 from eg, (44) and (45) 


(44) V- = 


(45) 


i/« 
’ 4 


r e x i x 

L /r~ 


erf /x”x 


ft 

4 


1 

.M* 


e X 2 x 


erf 


/v] 


/I7 


, e X 2 x . 
erf /Lx + “ erf 




X 2 X J 


with * i+1 and i 2 - 1-1 , 


One point is worth mentioning concerning fcfe algebra of complex 
numbers. If we take the root of a quantity like 

(D7) i + 1 s /2 e 4 ', i-1 = /2 e 4 . 


*) sold as "TRS80 Microcomputer" by Radio Shack in the U.S..This 
is not a commercial for microcomputers, indeed, the memory and 
speed of this computer are not very satisfactory. In the mean- 
time, however, more powerful microcomputers are available. 
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Program Unking of 
tho aubvou linos 


"SOim 1 " 


1 0 1* "SUtl?l" 
10*lH'271«.0?0P2 
S0764 

lost Ht 281 *. 04226 

P0 1 * 4 S 

1.0<ltH>si'*,0092'? 

05272 

lOStM'SCU*. 00015 
2014S 

lUfe.'n^sn*. 00027 

107iH-3^«,0000-! 

50&S8 

IGEbE-a+AtSTSiR 

+H<S8HR*R+ft 

»'29mR a $HA 

(SOj^cR-^nn 

aiS+tR'SMA 

"•16 

109tE“l-E 
110tN*ltH»Qils«0 
111'^SIH 2IRtC* 

COS SIR 

UStT»4RPtU»»2ftC« 

MaJRS 

list MeEXP CflfHB 
“I IN 

Jldsgc.B+Ti.tEMR 
‘B ‘I ’ 

1 IStPai L“P+U4+*1 

U • r. 

i ltsiC'n* in , 5 *mh I'l 
“1 N" 2 
1 17IH- H+Rt K*l •♦•0 » 
IP 11+1 

nt* ir hpg p.-je»- 

•lv..TC> MS 
11 V* IP HfS C'Mt- 
•IbOTO M S 
11‘Oi T kENF <-U>! 
t *EP CRR1SM 

*gnp 

ismh=t+'1»ce+>;c- 

1 • ll+SH It i 

ll'/tL&T+i ISE+S H 

+ 2 i H. 

l.StPCTUPH 


"sunR 2" 

CUM '‘ HIP?" 
fOil"*Rh+ lltlm 
MTU «I P> 

PH? i II- R'OOOTO 2 

Ob 

:mnr j diet d* 

n*r 

io’SUf MOLET JJ° 

»-» 

.■HKiMuMpBiiii »o: 

l »o 

207 !P»r-+*. l g«-IM :'’£5 
t|l»‘2M+l h*Ii 
SO'I' t F*F'+ COS II 

?09ti<*Pi$I!l h 

21 OH *MFiL«L+Gt 
ti*IN 1 

21 M IF li 70010 a 
07 

212t M«EflP CRP-11 
>tT»Q*r)MW 
'.no.' 

2l3t OH+COS 2RM 
R T*l H 

21-ltBM'li‘nil 2l?r+ 
1 Ml. H 
JistPCTi.mi 
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Pig. 25 roots of complex numbers. 

wo take only the firr.i. of two possible roots (Sec Piq^.25) 


if. 

/TFT * VI e '• » V2 (cos ir/8 + i sin n/8) 
i^i 

(D8) i / I r T « VI o ; « "/2 (cos 3n/8 + 4. sin 3n/8)u 

» V2 (sinn/8 + i cos n/8) 

we proceed in the same way with the quantity -. 

With these relations and the subroutines, the evaluation of egs. 
(44) and (45) is a straightforward procedure. Wo end up with the 
real and imaginary parts of i*o., v lR and v n and with simi- 
lar quantities for v ? > i.o., V SR and v 0I . 

In addition to that, the computer calculates also the difference 
between v and its asymptotic form, the amplified Helmholtz insta- 
bility wave 

(09) v 1D - v t - v 1H ; v 2D « v 2 - v 2H 


ll 

H 

: i 

i; 

i i 


with 


1 H 


i/n 

~r 




l V 


‘'ll 


2 II 


4 


e 1 
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As before, wo will calculate the real and imaginary parts also 
of the quantities in cq. (D9) , i.e., v^ DR & v.| D1 ns we * 1 08 

V 2DR & V 2D1‘ 

Finally, the program calculates the modulus of v^ R , which is 
equal to l v 2ii ^ nnd ^ha 1 2 3 4 "°^ u l us of v^ and v 2 , respectively . Since 
the memory of the microcomputer is quite limited, the mean pro- 
gram to calculate all those quantities cannot contain both sub- 
routines. Wo have two programs, "v-PRO 1 " with the subroutine 
"SUBR 1" for small and middle values of X and a second, almost 
identical program "v-PRO 2" with the subroutine "SUBR 2" for 
largo values of x. It is suggested to switch from "v-PRO 1" to 
"v-PRO 2" at values of X between 5 and 8. it has been chocked, 
that also the series expansion (D1) produces overlapping values 
for the velocities up to x * 2*3. The operation of the programs 
starts with the computer asking for 

(i) the initial values of x ("x- START?") 

(i,i.) the stepwidth Ax for the consecutive calculation of a 
series of x values ("SCHRITT?") 

(iri) the final value of x of the data set ("x-ENDE?") 

The programs "v-PRO 1" and "v-PRO 2" are listed on the next pages. 

There is another version of these programs which provides data on 
casette tape for further use of the following programs calculating 
u^ {Appendix E) . Since only v 2QR and v 2DJ are the necessary quan- 
tities to bo supplied for these further u-programs, the v-program 
is much shorter, so that both subroutines for small and middle x 
and for large x can be comprised in the program. Each data set 
stored on the data tape consists of 4 words, i.e. 

1) the label "v 2D " 

2) the value of x 

3) the value of v 2DR 

4) the value of v 2D .j 
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Program listing 
"v-PRO 1" 


2s INPUT "X-STA 
RT?"» A<45>* " 
SCHRITT?"»A< 
46>i "X-ENDE? 
"»A<47> 
3sX=A<45> 

9i IF X>OGOTO 1 
1 

10iA<37>«0ift<38 

>=0iA<39)*3» 

A<40)=08GC»T0 

27 

ll80=T<X*42>lU= 
COS <R/8)iV= 
SIN W8) 
12iR=0U«I«0V 
I3JGC5SUE 101 
15sA<33>=AU+BV 
16sAC34>=-AV+BU 
178R=OV:I=OU 
18:G0SU3 101 
20:A<35>=AV+BU 
21 sA<36)=~AU+BV 
22s Y=-<-TH>/'<4*7 
<T2>> 

23sA<37>=Y*<A<3 

3>+A<35>> 

248 A<38)=Y*<A<3 
4>+A<36>> 

258 A<39>=Y*<A<3 
4>-A<36>> 
26SA<40)=Y*(-A< 
33>+A<35>> 
27!Y=-<TI0*<EXP 
x>/< 4 *«ra 2 » 
288AC41)=A<37>- 
Y+COS <X-<R/ 
8 ?) 

29sA<42>=A<38>- 

y»sin <y,-a/ 

8 >> 

30sA<43>=A<39>+ 
Y*COS <X+<3H 
/8 >> 

31sA<44>=A<40>+ 
Y+SIH <X+<3J1 
/ 8 » 

328 PRINT "X="! 
PRINT X 

33s PRINT "V2R=" 
SPRINT A<37> 
348 PRINT "V2I«" 
SPRINT AC38) 
35 8 PR I NT "V1R=" 
SPRINT A<39> 
36 8 PRINT "VI I=" 
SPRINT A<40> 
37s PRINT "V2DR= 
"SPRINT AC41 
> 

38s PRINT "V2DI= 
"SPRINT A<42 

> 

398 PRINT "V1DR= 
"SPRINT A <43 
> 


PRINT -Y 
42SPRINT "ABS<V 
2>="sPRINT f 
<A<37>+A<37> 
♦A<33>*A<38> 

> 

438PRINT "ABS<V 
1>="8PRINT f 
<A<39>*A<39> 
+A<40>*A<40> 

> 

44sX=X+A<46> s IF 
X<=A<4?>GflT0 
9 

45 s END 
101!' B SUBR1" 
102sA<27)=. 07052 
30784 

103sA<28)=. 04228 
20123 

104SA<29)=, 00927 
05272 

1058A<30>=. 00015 
20143 

1065 AC31 >=.00027 
65672 

10?:A<32>=. 00004 
30638 

108sE=<l+A<27>*R 
+A<28>*R*R+A 
<29> 1 * ! <R' s 3>+A 
<30>*<R a 4.HA 
<31 )'KR A 5>+A 
<32>*<R A 6) > A 
~16 

109SE=1-E 
110814=1 ^H=08K=0 
Ills S=SIN 2IR8C= 
COS 2IR 

112sT=4RR8L=2RC8 

M=2RS 

1 13! W=EXP <NI>8D 
=NN 

ri 4 ! 2 =ai+T)*EXP 

<H/ 4 > 

115!F-<L-R*<W+1/ 

W>>/2 

U6ii?,=<M+.5*N*<W 
-iA‘\y>'Z 
1 17:H=H+PsK=K+Q! 
N=N+1 

118! IF ABS P>1E- 
4G0T0 113 
1198 IF ABS Q>1E- 
4 GOTO 113 
12Q! T=EXP <-II>: 
L=EXP <RR>8M 
=2J1R 

121IR=T*<LCE+<C- 
1 VN+2H/I> 
122! B=T+<LSE+S/H 

+2! jo 

123: RETURN 
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It "V-PROl" 


408 PRINT "VI Bis 
"SPRINT A<44 

> 

41 8PR1NT "ABS<V 
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Program listing 
"v-PRO 2" 


li "V-PROS" 

21 INPUT "X-STA 
RT?"* A<45>» " 
SCHRITT?">A< 
46>» "X-’ENBE? 
"fA<4?> 
3lX*A<45> 

91 IF X>OGOTO I 
1 

10iAC37)=0iA<38 

>=0iA<39>>=0! 

A<40)=OlGOTO 

27 

1II0=T<W2>:U= 

cos <n/8>*v= 
sin <n/e> 

12IR=0U: I=0V 
14IG0SUB 201 
13sA<33>=AU+BV 
16!A<34>s-AV+BU 
171R=0VH=0U 
19IG0SUB 201 
20«A<35>=AV+BU 
2HA<36>=-AU+BV 
22tY«-an>/'<4*7 
<F 2 >> 

23IA<37>«Y+<AC3 

3>+A<35>) 

24iA<38)=Y*<A<3 

4>+A<36>)> 

25sAC39>=Y+<A<3 

4>~A''36)) 

26! A<40)=Y>K-A( 
33>+AC3S>> 
27* Y--<T JO*<EXP 
X>/<4*fa2>> 
28iAC41?=A<37>- 
Y*CCJS <X-<V 
8>> 

29iAC42)=A<38)- 
Y*SIN <X-<R/ 
8>) 

30lAC43>=A<39>+ 
Y#COS <X+(3J! 
/8» 

3l!AC44)=A<:40)+ 
YiS IN <X+<3JI 

321 PRINT "X="! 
PRINT X 

33! PRINT "V2R=" 
! PRINT A<37> 
34s PRINT "V2I-" 
SPRINT A<38> 
35s PRINT "V1P=" 
.‘PRINT A<39> 
38s PRINT "VI is" 
SPRINT A 1 40? 


371 PR I NT "V2DR= 
"« PRINT A<41 
> 

38SPRINT "V2BI= 
"SPRINT A<42 

) 

39! PRINT "VIBRs 
"SPRINT A<43 

> 

40* PRINT "V1BI= 
"SPRINT A<44 

> 

41SPRINT "ABS<V 
2H=VlH>="s 
PRINT -Y 
42SPRIHT U ABS<V 
2)="!PRINT «T 
<A<37)fA<37) 
+A<38>*AC28> 
? 

43* PRINT ”AB3<V 
1>=" SPRINT ■! 
<A(39?*Af39) 
+A<40>*A<40) 
> 

44*X=X+A<46>! IF 
X<“AC47>G0T0 
9 

45:ENIi 
201 s "8UBR2" 

202s G=RR+I I : H= 
ATN CI/R> 
203* IF R>OGOTO 2 
06 

204: IF I >OLET B= 

n+n 

205! IF KOLET D= 
B-Jt 

206!M=l!P=l!K'=0! 

L=0 

207 : F'=F'+ f , . 5-M > .'Q 
sN=?2N+l >=PB 
208!F=P*C0S N 
209sG=P*SIN N 
210!K=K+F!L=L+G! 
( 1 = 11+1 

21 is IF M<7GOTO 2 
07 

212:M=EXP CRR-II 
):T=Q*T7t:H=T 
<JIQ) 

21 3* A=N*COS 2RI- 
R/T-K.‘H 

214!B=M*SIN 2RI+ 
I/T+L'H 
21 5! RETURN 
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Besides that output feature, the "tape output version "v 2D “PRO" I > 

is identical to v-l’RO 1 4* v-PRO 2f the switching from the j 

expansion for small and middlo x (SUDR 1) to largo & (SUBR 2) I \ 

takes place at S ■ 6. Finally, wo provide in Table 2 a data j 

survey for the different v-velocity constituents. I 

I « 

1 1 

I f -* 

* * . 

k 4 
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Program linking 
Hv 2D^ R0 " 
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1 1 41 V2D-PK0" 

2« INPUT •‘X-STA 
RT?%A<2?>»* 

$chritt?%a< 

28>*"X-ENDE? 

"iA<29> 

3*N«A<2?> 

9* IF X5QG0T0 i 
1 

lO*Ar,X05-0«A<3t 

5»0«A<32>»0* 

A<3S5»0*GOTO 

2*j» 

itiO«f<X*T2>*U« 

eos <*ve>«Y« 

SIN Ok‘05 
12* R»0U* I*0V 
13* IF X< 6G0SUB 
101 

14* IF XW'GOSUB 
201 

l5*A*;34i*AU*IV 
I6«A«5)*“AV4*U 
lTiR«OVU«OU 
IQ* IF X<6G0SUB 
101 

HI IF X>»iG0Sl)B 
201 

SOiACMJ-AVFBU 

21UH37)«-AU+BV 

22*Y«-ajO/<4*l 

<«» 

23»A<30>«YHA<3 
4>+A<36>> 
24tfK31)*¥*<A<3 
3>4.WS?» 
23*fH32)*Y<KA<3 
35-A<3?» 
26*fH33)«Y*<-rt< 
34)*fK36>5 
S?«V»-em«EXP 
XV<4*W25> 
30tACS«5*AC305** 
YiCOS Oi-CR' 
8?) 

3HA<57)-ACSt)- 
YtSIN <X-<R/ 

m 

32* A<555“X 
33«PR1NT i*V2D* 
IA«S> 

44«X«X*A<28)*1F 

X<-A(295GOTO 

9 

45* END 
101*"SUBR1" 

1 02 * K“ . 070323076 

4 

103*L«. 042262012 
X 

104* M*. 00927052? 

2 

l 05* N“. 000132014 
X 


106* P«, 0002T636? 
2 

10? IQ-. 00004X08X 
8 

109iE*«1*KR*LRR4 
M+*R'3HN*tR 
*45+P*<R'3>4 
G+«R*6> 5**16 
109*E»1-E 
110iN»liH"0*K“0 
lll*$*SiN 2IR*C- 
CGS 21R 

1 l2*T“4RRlL»2RC* 
M*2RS 

113iW»EXP <Nl>lD 
*HN 

il4»2*‘«ll+T>*EXP 
CD 45 

113*P«vU-R+CM*l^ 

W*V2 

H6*«»*m.3+N*<W 

ll?iH*H+P»R*K40t 

H-NU 

US* IF ADS PMt- 
4C-0T0 113 
1 19* IF ADS CSME- 
4 KOTO US 
ISO* T»EXF <~I IV* 
L»EXP CRR>*H 
*2HR 

121iA»t*aCE+CO 
n N*2ivn5 

122* D«TiU,SE*S'H 
+2WJ15 
123* RETURN 
201 t “SUDR2" 

202* Q-RR* Mill* 
ATN aH«> 
SOS* IF R50G0T0 2 
06 

204* IF IMH.ET D» 
D+Jl 

203* IF HOLEY D« 
D“ll 

206*t'M*P*l*k"Q* 

L«0 

207*P«P+C.5-N>^Q 

iN*v2M*-t5*D 

2Q8*F**P*C0S N 
2G9*»P+51N N 
210*K-R4F«l.»UG* 
N»N+l 

21H1F CK7G0TO 2 
0? 

212* M“EXP CRSMI 
)*T-«Tf.*UH«7 
0105 

21X*A*tWCIS 2RI“ 
R-’T~k»H 

2I4*B«H*SIN 2RI+ 
l/T+L^H 
213* RETURN 
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u~computations at y < o, tables 

The basic ideas of the u 2 field computations have already been 
outlined in section 3, The calculation proceeds as follows (see 
Fig. 26 ) 



We calculate the velocity u 2D , induced by the v 2d velocity distri 
bution. u 2D is computed for a set of, say, 15 field points with 
varying -y and at a constant x location. We are summarizing the 
induced velocities from each (source distribution) interval 65 , 

We start the computation at sufficiently large negative values 
of 5 to have a sufficiently low value of v 2D * We use an integra- 
tion procedure utilizing trapezoidal (source) elements between, 
say, 5 = -6 and 5 = o. For the real part of v 2D , which has an 
infinite slope at { = c (see Fig. 6 in section 2) , we use for 
the v 2R contribution to v 2DR a parabolic piece of curve between 
5=o and 5 = o. 1 . For this region we use the analytic solution 
eq. (82). For the region between 5 = o.l (or % = o for the 
imaginary part) and 5=6, we use again the trapezoidal element 
model. The stepwidth for our calculation is 65 = o.l. Above 5 = 6, 
the real part is set equal to zero (no contribution above 5 = 6) 
and the imaginary part, which has a -1 /4 /g behaviour, is evalu-^ 
ated with the analytical far field solution, eq. (79) . The 


computation of u 2BK and u 21)V iti done simultaneously. The valueu 
of v 21) for A • o are computed by Lho u 2 - prog ram, bocauao it 
consists just of the (negative) Holmholta solution, Pol* values 
between A c o and A » 6 the values aro taken from a data tape 
produced by the program discussed in Appendix p. For values 
above ■*, « 6 we use again asymptotic solutions, ns mentioned 
above. 

Since many expressions in equation (79) for the far field and 
In equation (02) for the field at the plate edge are the same, 
we use the sawn subroutine , i*e,, instructions 2011209 of the 
program, 'the trnpeaoldal integration segment is the subroutine 
between Instructions 101 and 105. 'l'he field calculation between 
i. a o and t “ o,1 is done, between instructions 60 and 65, and 
the far field calculation is done between instructions 50 and 52. 
'l'he rest of the. program consists more or less of program manage- 
ment like loop controls etc, as well as the calculation of the 
llelmhoita wave field which has to bo added to the values of « 2D 
to obtain u 2 * m addition, the reading of the v 2)3 values from the 
data tape is cheeked with the p. of each data sot read from the 
tape against a /, generated internally by the program. For a 
given x of the field points, the program calculates 15 values of 
-y with a step width of 0.05, 'this takes about 6 hours (I), so 
the data set for u-, consisting of 16 x-mibsobs took about 96 
hours. This is outrageous, but the computing time is free for a 
microcomputer and much time As saved during the program estab- 
lishment phase, since the microcomputer can be used continuously. 
The BASIC computer program for u ? is given on the next page and 
it is followed by a set of Ug field data (see Table 3 ) . Plots of 
those data are given in Pig. 10, 
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"u 2 -PRON" 
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lfya-PROH" 

2: INPUT 
X 

3:A*-6.1:N«2?: 
11*46! V=l.E-8 
:E=0:F*0:3=< 

) 

4! A=A+. 1 
5s IF AHOGOTO 
10 

6:W=2+EXP A:G* 
W*C0S fA-OU* 
8»:H*M*SIN 
cA-a^eio 
7 i GOTO 14 
10: INPUT S"V2B" 
1 A>: 65 ) 

li:G*A(I66>:H*A< 
67 i 

12: IF fl<>A<65> 
PRINT "LESEF 
EHUER" 

13! IF A*. 1G0TO 
60 

14:10= <G-EH10:L 
*E 

15SG0SUS 101 
16:ACN>i»A<NV+U 
l?:K-i.H'“F?»lO:L 
=F 

18! GOSUB 101 
19«AiM)»ACM>+U 
20: IF A*6. OGOTO 
50 

21:Y“Y+.05:N=N+ 
1 : 11 = 11+1 

22: IF Y<0. 8G0T8 
13 

2?:E=G;F=H: Y=0 
24 :H=2T: 11=46 
25: IF fl<‘ 6G0TD 4 
26: PRINT "X=VX 
2?: PRINT ‘'Y-">- 

Y 

28 !Us 2 '*e;:p cx-y 

2?!ps-HfN>+W* 

81 N ' l+Y-a 

§ ."i 

30: inWlHHN* 
COS >:'+Y-'.U 

s ■ > 

21 i PRINT "U2P=“ 

> R 

32! PRINT "1.121*“ 

.1 

23ipFINT "ABS'U 

2 '=", f.pp+Il 

34: IF RbOLET R= 

V 

35: C-hTN > i -rj 


36: IF R OGOTO 3 
■s 

37: IF rOLET 0= 
0+K 

30: IF KOLET 0* 

o~n 

39 SPRINT “PHASE 
>U2‘'£EG="> CO 
jp+iso 

40SY*Y+.05!N*N+ 
1 : 11 = 11+1 

41: IF V<.8G0TCi 
26 

42! END 

50: GOSUB 201 

sIwHhacmhci 
•k-o+csin <d 

52: GOTO 21 

60s G =2+ < EXP AH 
COS- ih-CK+ 8) 
':l *:G-EH10 
:L=E 

61 : GOSUB 101 

62!A<N''*Atl'D+U 

63: GOSUB 201 

64ifi+NT=A(NH<< 
•TWHUHCO+T* 
SIN U5-25HC 
<2+-r.l>/]piG 
sG-T. 1 

65 .’GOTO 17 
101 : £=h-. 1 ! Y=Y+V 
ICCsPaOl-BHCX-B 
.HYY 

103iC!*<.X-AH<X-A 

>+YY 

i045U=ii;"jn*cY+< 
ATN uA-XHY 
>-hTH U'B-X) 
•Y>‘i-i.iL.. J <K+ 
V>>+M~BH.5* 
LN fO'PH. 1> 
: Y*Y-V 
105: RETURN 
201;P*hTN O.VCX+ 

202: IF X OLET D= 
B-f it 

.'03:U=T> '” !*YV t+V 
2f4i';-2<TiHW'>M 
CO? ij 2' 
20?:! =2"* ■. f« AW'H 
SIN • B 2» 
506>.TsttTH I'l'ill- 

207: IF • M-h.( 0 
SET TaT+R 

» cos vn 

2 1 . • EH * : N+ 
H+C : 1 U+Hrrt* 
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